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. Introduction  and  Summary  of  Results 


The  problem  of  propagation  of  vaves  in  a random  medium  appears  in  a 
number  of  areas  of  research  and  applied  science.  Some  examples  are  atmospheric 
optics,  radio  astronomy  and  underwater  sound.  The  problem  is  furthermore  sm 
old  one  which  has  been  studied  extensively.  The  earlier  work  (summarized  in 
the  monographs  of  Tatarskii  and  Chernov  ) employed  the  Rytov  approximation. 

In  this  approximation  the  logarithm  of  the  amplitude  is  computed  using  first 
order  perturbation  theory.  The  Rytov  method  is  applicable  whenever  the 
Intensity  fluctuations  are  small.  When  the  wavelength  is  small  it  reduces 
to  first  order  geometric  optics  or  WKB.  More  recently,  a different  approxi- 
mation which  reduces  the  problem  to  a Markov  process  has  lead  to  considerable 

progress  in  cases  where  the  intensity  fluctuations  are  not  small.  This  method 

3 

is  explained  in  Tatarskii's  second  book  and  in  two  excellent  reviews  of  the 

5 

recent  literature.  ’ Nevertheless,  important  problems  remain.  In 
particular,  there  does  not  exist  a global  view  of  what  is  going  on  in  the 
so-called  saturated  regime  where  the  intensity  fluctuations  are  important. 

In  this  paper  Feynman's  path  integral^  is  applied  to  the  problem  of 
wave  propagation  in  a random  medium.  It  provides  a natural  and  systematic  method 
for  attacking  the  problem,  especially  when  the  intensity  fluctuations  are 
large  and  the  Rytov  approximation  fails.  The  path  integral  is  widely  used 
in  quantum  mechanics  and  statistical  mechanics  but  it  is  expected  that  many 
readers  will  not  be  familiar  with  it,  thus  the  paper  is  meant  to  be  self- 
contained.  The  reader  who  desires  further  background  information  on  path 
integrals  will  do  well  to  consult  the  book  of  Feynman  and  Hibbs.^ 
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Because  some  readers  wil^.  not  be  familiar  with  path  integrsJ.s  there  are  i 

some  pecularities  in  the  organization  of  this  paper.  In  real  situations,  i 
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frequently  have  a power  law  spectrum  in  the  scale  size  of  fluctuations. 

Path  integrals  are  capable  of  handling  all  these  complications,  (in  fact 
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the  author  first  developed  the  method  for  propagation  of  so\md  in  the  ocean  , 

a problem  vdiich  has  these  complications  and  more.)  However,  it  is  veistly 

easier  to  explain  the  path  integral  method  for  an  idealized  medivan  which  is 

statistically  homogeneous  and  isotropic  and  whose  fluctuations  eire  charac- 

terized  by  a single^  typical  scale  size  L ( small  compared  to  the  distance  R of 

propagation).  The  bulk  of  the  paper  is  therefore  devoted  to  a study  of  this 

idealized  situation.  Once  this  has  been  done  the  transition  to  realistic 

media  is  relatively  simple.  However,  this  manner  of  presentation  has  a defect 

for  which  only  an  apology  can  be  offered.  Because  of  the  temporary  restriction  ' 

I 

to  a single  scale  size  L,  results  which  are  directly  applicable  to  atmospheric 
optics  do  not  appear  until  late  in  the  paper  (specifically.  Secs.  (7)  and  (8)). 

Finally,  to  illustrate  the  power  of  the  path  integral  method  (and,  hopefully,  I 

motivate  the  reader),  a number  of  results  for  the  idealized  problem 

will  be  summarized  below.  The  translation  of  these  results  to  more  complicated 

cases  is  generally  straightforward:  the  details  are  given  in  the  text. 

Listing  the  results  will  require  the  definition  of  some  symbols.  This  will 

(temporarily)  be  done  in  terms  of  the  idealized  problem  and  the  reader  who  j 

I 

has  worked  on  propagation  in  a random  medium  will  find  that  they  are  familiar  i 

objects;  e.g.,  Tatarskii's  phase  structure  function  D.  For  other  readers, 
the  motivation  for  these  definitions  will  become  apparent  in  Secs.  (2)  and  (3).  I 

Actually,  there  are  two  distinct  kinds  of  problems  of  propagation  in  a | 

random  medium,  corresponding  to  whether  the  scattering  angles,  single  and/or  j 

multiple,  are  large  or  small.  If  the  fluctuations  are  weak  so  that  a single  i 

scattering  approximation  (Born  approximation)  applies  there  is  little  dis- 
tinction between  the  two  cases.  However  in  a multiple  scattering  regime. 
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wlilch  Is  the  case  of  Interest  here,  the  two  kinds  of  problems  are  very- 
different.  This  Is  Illustrated  in  Fig.  (l).  The  considerations  of  this 
paper  will  be  restricted  to  situations  i^ere  the  single  and  multiple 
scattering  angles  are  small.  This  is  sufficient  to  cover  the  applications 
mentioned  above.  The  large  emgle  multiple  scattering  situation  is  like  a 
problem  in  radiative  transport  and  is  most  efficiently  treated  by  other 
methods . 

It  will  be  assumed  that  the  problem  ceui  be  reduced  to  a scalar  wave 
equation  with  an  index  of  refraction  n(x,t)  which  may  depend  on  the  frequency 
to  = ck.  In  a homogeneous  medium  { n ) is  a constant  and  for  waves  of  a fixed 
frequency  can  be  set  equal  to  unity.  Defining 

p(x,t)  = 1 - n(x,t)  (1.1) 

W will  be  taken  to  have  a zero  mean  and  a covariance 


(^l(x,t)  u(x',t')>  = p(lx  - x'l.t  - t'). 


(1.2) 


Q 

It  will  be  further  assumed  that  either  u is  a Gaussian"^  random  field  or 
, / 2,  i 

that  kLiu  ) is  small,  in  which  case  the  distribution  need  not  be 
specified. 


Let  the  two  dimensional  vector  r = (x  , y ) label  the  location  of 
30  o o o 

a point  source  in  the  plane  z = 0.  Then  in  a plane  of  constant  z > 0, 

the  signal  will  be  E(z,r,r  ,t)  where  r s (x,y)  specifies  the  transverse 

o 

coordinates  of  the  observation  point.  The  total  range  of  propagation 
will  be  denoted  by  R and  for  |r[ , |r^|  « R and  a CW  source  it  is  useful 
to  define  a complex  envelope  S by 


E(z,r,r  , t)  = Re 
o 


f -»  -»  i 

(S(z,r,r^,t)  e 


(kz  - u)t) 


(1.3) 
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The  time  dependence  of  (1  is  due  to  fluctuations  in  the  medium.  It  will 

be  assumed  that  the  full  wave  equation  for  £ can  be  approximated  by  the 

3-7  11 

parabolic  wave  equation  ’ 


(‘i 


^ 7^  - ku(r,z.t)j5(r,r^,z,  t)  =0 


(1.4) 


.2  ,2 

O ^ pi 

where  V = + plus  a boundary  condition  at  z = 0 


-♦  (4kz)  ^ exp 


ik(r  - r ) 
o 


2z 


(1.5) 


If  L and  T are  the  characteristic  lengths  and  times  over  which  \x  changes, 

the  validity  condition  for  Eq.  (1.4)  are  (i)  kL  » 1,  (ii)  kL  « u)T  and 

2 ) 

(iii)  that  the  rms  multiple  scattering  angle  ((^  )R/L)  should  be  small. 

Feynman's  path  integral  gives  the  solution  to  the  parabolic  wave 

equation  in  terms  of  a (strictly  speaking)  infinite  dimensional  integral. 

It  turns  out  that  this  integral  can  be  studied  in  almost  exactly  the  same 
12 

way  as  Mercier  originally  attacked  the  phase  screen  integral.  The  result 
is  that  propagation  in  a statistically  homogeneous  medixom  is  very  similar 
to  the  phase  screen  problem.  This  will  continue  to  be  true  in  rather 
general  inhomogeneous  media,  of  which  the  phase  screen  is  a special  case. 

In  order  to  indicate  what  can  be  learned  from  the  path  integral  it  is 
necessary  to  review  some  known  features  of  propagation  in  a random  medium. 
The  qualitative  character  of  £ is  determined  by  two  parameters  * and  Q 
defined  by 

R 


= k^<^|^^ti(e^  z , t)dzj 
00 

= k^Ry’p(|zl,0)dz  + 0(L/R) 


(1.6) 
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where  e Is  a unit  vector  in  the  z direction  and  in  the  second  line  it 


z 

has  been  assumed  that  R » L and 


2 

n = 6kL  /R 


(1.7) 


The  parameter  § is  just  the  rms  phase  fluctuation  as  computed  in  first 

1 

order  geometric  optics  and  serves  as  a measure  of  the  strength  of  the 
fluctuations.  The  other  parameter  Q is  essentially  the  square  of  the 
ratio  of  the  scale  size  L to  the  extent  of  a Fresnel  zone.  As  shown  in 
Figure  2,  if  ^ is  less  than  one  or  less  than  n,  then  the  Rytov^  ^ 
approximation  is  valid.  In  the  region  where  the  Rytov  approximation  is 


valid,  the  problem  can  be  considered  to  have  been  solved  years  ago. 

The  intensity  fluctuations  (scintillations)  are  small  and  the  relation 
between  S and  u is  simple  and  direct.  Also,  as  shown  in  Figure  2,  when 
both  4 and  4 /Cl  are  greater  than  unity,  the  fluctuations  in  S saturate. 

In  particular,  the  variance  of  In  |£|  = In  I approaches  a constant  of 

order  unity  and  the  properties  of  S are  determined  more  by  statistical 
considerations  than  by  the  detailed  properties  of  p. 

Path  integral  methods  have  nothing  new  to  add  when  the  Rytov  approximation 

is  valid.  The  considerations  of  this  paper  will  therefore  be  restricted  to 

the  saturated  regions.  There  is  then  a small  parameter  ot  = £1/#  whose  order 
lU 

of  magnitude  is 


or 


6L 


3/2 


(1.8) 


The  path  integral  allows  the  calculation  of  any  moment  of  S as  an 
asymtotic  series  in  a.  The  result  is  that  S is  uniformly  distributed 
in  phase  and  that  the  monents  of  intensity  I = ((J|  are  given  by 

<l")  = n;<I>"[l  + in(n  - l)Ckr  + 0(a^)]  (1.9) 

where  C is  a calculable  constant  of  order  unity  whose  precise  value 
depends  on  the  spectrum  of  (i.  In  the  limit  or  = 0 the  distribution  is 
therefore  Raylelgh^^  with 
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P(I)  = ^ .xp[-  ^ 


(1.10) 


However,  the  correction  grows  with  n and  cannot  be  neglected  for 
n > (2/0)^.  It  follows  that  there  must  be  significant  deviations  from 
a Rayleigh  distribution  when  I/(l)  is  greater  than  ~ (2/a)^. 

In  addition  to  the  distribution  of  S,  one  also  wants  to  know  the 

coherences  in  space  and  time.  Recent  work  on  coherences  has  been  greatly 

3-5 

facilitated  by  the  observation  that  under  certain  conditions  the  prob- 
lem can  be  replaced  by  a simpler  local  or  Markov  one  where,  in  effect, 
one  makes  the  replacement 


with  r = (x,y)  and 


p(x,  t)  -♦  6 (z)p  (1  r]  , t) 


p(fr(,t)  = + z^)^,tjdz 


(1.11) 


(1.12) 


2 2 ^ 

Note  that  ^ is  equal  to  k Rp(0,0)  and  it  will  be  convenient  to 
use  the  function  3 to  define  T aad  L by  the  expansion 


„/  -»2  2 

2 . *2/  r t 

k Rp(|r|,t)  = * 11  - — + ... 


„ 2 „ 2 
2L  2T 


(1.13) 


Within  factors  of  order  unity,  the  L and  T so  defined  will  be  equal  to 
the  length  and  time  over  which  the  original  covariance  p is  non-vanishing, 

3-5 

It  has  been  pointed  out  by  several  authors  that  in  the  Markov 
approximation  the  coherence  of  0 and  S can  be  computed  exactly.  It  is 


((S  (r ',r  ,t')^(r,r  , t)> 
o o 

S\r',T')<S  (r,r  ) 


= exp  [-iD(r  - r',r  - r',t  - t')]  (1.14) 

o o 


where 


(S  (r.r  ) = <I>^  exp  ik^^— ~ 

O O 2R 


(1.15) 


and  D Is  the  phase  structure  function  of  first  order  geometric  optics 

1 

D(r,r^,  t)  = 2k^R^[p (0,0)  -p(lur+  (1  -u)r^|,t)]du  (1.16) 


The  phase  structure  function  always  appears  in  an  exponential  and  in  the 
saturated  region  where  $ is  large,  D can  be  approximated  by  an  expansion 


in  r,r  and  t 
o 


-♦2  -»2  ->  -» 
„r+r+r"r  2 

D(r,r,t)=$  ^ 

3L  T 


(1.17) 


Coherences  are  then  characterized  by  two  parameters  $/T  and  i /L, 

The  literature  is  somewhat  confusing  as  to  the  validity  conditions 

for  Eq.  (1.14).  It  turns  out  that  the  approximation  leading  to  Eq.  (1.14) 

has  a very  simple  interpretation  in  the  path  integral  formalism.  In  the 

next  section  is  will  become  evident  that  for  the  isotropic  medium  under 

consideration  Eq.  (1.14)  is  valid  as  long  as  the  parabolic  wave  equation 

is  valid.  From  the  path  integral  one  can  actually  compute  the  first 

correction  to  Eq.  (1.14).  It  is  of  order  of  the  rms  multiple  scattering 

i 2 J 

angle  (R/L)^(u  ) which  must  be  small  if  the  parabolic  wave  equation  is 
val  id . 


For  small  a the  path  Integral  also  allows  the  calculation  of 

(<J*(r,r  , t)(S’^(r r ',  t ')(f  (r  ' r ' t ' ')(S(r  " r ^ t ' ' ^))  and  more  generally 
o o o o 

an  arbitrary  2n-th  order  moment.  In  the  limit  or  = 0,  the  real  and  imaginary 
parts  of  (f!*  are  Jointly  Gaussian.  To  see  the  use  of  this  result,  let  us 
consider  a typical  question  of  practical  interest.  Take  a fixed  source  and 
receiver  so  that  S'  is  a function  only  of  time  and  suppose  that  at  t = 0 


iT/d'  is  known  to  have  a value  n-  An  interesting  practical  question  is 
o 

then  what  is  the  probability  P(n')  that  C.{t)/C^  will  take  on  the  value 
n'.  Since  ^ has  a Gaussian  distribution,  P(n')  is  simply 


P(Tl') 


exp 


V- 


1 - e 


-D(t) 


:(l  - =-“<•>) 


(1.18) 


where  D(t)  = D(0,0,t)  ^ (t/T)  . The  qualitative  behavior  of  is 

indicated  in  Figure  3.  It  is  evident  that  the  signal  stays  in  one  quadrant 

of  the  complex  plane  and  is  therefore  coherent  over  a time  of  order  T/^ . 

A further  property  of  Gaussian  statistics  and  a covariance  of  the  form 
2 

exp  [-^('^t/T)  ] is  that  the  signal  will  move  in  a straight  line  for 

times  less  than  ~ T/$ . One  can  ask  the  more  general  question  of  given 

that  5(r,r  ,t)/o  (r,r  ) is  equal  to  what  is  the  probability  that 
o o o 

<S(r',r^,  t')/{S  (r',r')  will  be  equal  to  Tl',  llie  result  is  just  Eq,  (1,18) 
o o o 

with  D(t)  replaced  by  D(r  -r,r  -r,t-t). 

o o 

As  stated  above  the  Gaussian  statistics  leading  to  P(ti')  are  obtained 
by  computing  moments.  Again  the  approximation  scheme  breaks  down  for 
moments  of  order  (2/a)  ' and  Eq.  (l.l8)  is  valid  only  for  |ril  and  |n'i  less 
than  (2/a)^^**.  Actually  the  order  a corrections  to  any  moment  are  calculable. 


They  are  most  important  for  intensity  correlations  where  they  lead  to 
coherence  tails  of  order  a which  are  small  but  fall  much  less  rapidly  than  e~^. 

The  path  integral  also  provides  a simple  method  for  calculating  the 
correlation  between  waves  of  different  frequencies.  In  the  saturated  region 
where  $ > 1 the  result  is,  for  |a)-(i)'l  small  compared  to  uT  = ^u+{o') 


K&  (a)' ) (u))^ 


(1.19) 


8 


> 


» 


» 


» 


where  is 
g 


R 


o 


and  for  a single  scale  medium 


(1.20) 


(1.21) 


—2  “2  ^ 

with  u = c R p(0,0)  where  c is  the  unperturbed  group  velocity.  For  a non- 
o 2 S 

dispersive  medium  When  a is  very  smell  the  second  factor  on  the 

right-hand  side  of  Eq.  {I.I9)  falls  much  more  rapidly  than  the  first  one.  The 


first  factor  exp 


. .2-1 

- - 1 

3 

1 

3 

2 1 

1 0)  / 

^ g ' J 

can  then  be  replaced  by  unity,  in  the 


limit  a = 0 the  higher  order  correlations  in  frequency  are  Gaussian.  One  can 
then  obtain  probability  distributions  in  frequency  from  Eq.  (I.I8)  with  exp  [-D/2] 
replaced  by  the  right  hand  side  of  Eq.  (I.I9)  eind  e~^  replaced  by  its  absolute 
value  squared.  It  is  worth  noting  that  a first  order  geometric  optics 
calculation  misses  the  second  and  dominant  factor  on  the  right  hand  side  of 
Eq.  (1.19)  and  therefore  vastly  overestimates  the  range  of  coherence  in 
frequency. 

It  can  be  seen  from  the  path  integral  that  saturation  corresponds  to  the 
appearance  of  multiple  Fermat  paths  which  satisfy  a perturbed  ray  equation. 

The  signal  tends  to  propagate  along  these  Fermat  paths  and  because  there  are 
many  of  them,  they  interfere  and  produce  Gaussian  statistics.  They  will 
become  manifest  in  an  experiment  with  a pulsed  source  where  the  received 
signal  will  tend  to  show  several  arrivals.  These  multiple  Fermat  paths  are 

responsible  for  the  factor  A in  Eq.  (1.19). 
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With  one  exception  these  results  can  easily  he  extended  to  statistically 
inhomogeneous  or  anisotropic  media  and  to  media  with  multiple  scales.  The 
exception  is  that  A(uj)  defined  in  Eq.  (I.I9)  cannot  be  computed  for  certain 
multiple  scale  media.  Actually,  the  path  integral  yields  further  Information 

in  the  case  of  multiple  scale  media.  It  appears  to  be  only  partially  under- 

3-5 

stood  that  in  this  case  there  are  two  distinct  saturated  regimes.  An 
examination  of  the  path  integral  shows  that  there  are  indeed  two,  one  of 
which  (the  fully  saturated  regime)  is  analogous  to  the  saturated  regime  in 
single  scale  media  and  another  one  (the  partially  saturated  regime)  is  new. 

Many  experiments  in  atmospheric  optics  lie  in  the  partially  saturated  regime 
and  this  case  is  treated  in  some  detail  (Sec.  8).  The  fundamental  distinction 
between  the  fully  and  partially  saturated  regimes  shows  up  in  correlations 
between  waves  of  different  frequency.  In  the  fully  saturated  regime  the  real 
and  imaginary  parts  of  (w)  are  Jointly  Gaussian  random  variables.  For 
partial  saturation  (to)  acts  like  a random  phase  times  a Gaussian  object. 

A consequence  is  that  propagation  of  narrow  pulses  is  qualitatively  different 
in  the  two  regimes.  Depending  on  the  medium  there  may  be  further  qualitative 
differences  between  full  and  partial  saturation. 

The  detailed  organization  of  the  paper  is  as  follows.  Secs.  (2)-(6)  and 
Apps.  (a)  and  (B)  are  devoted  to  the  idealized  homogeneous,  isotropic  medium 
with  a single  scale  size.  In  Sec.  (2)  the  path  integral  is  introduced  and 
applied  to  the  calculation  of  the  first  and  second  moments.  App.  (A)  contains 
the  calculation  of  the  error  in  Eq.  (I.1I4).  Sec.  (3)  is  de  voted  to  the 
calculation  of  higher  moments  when  a is  small  and  Sec.  (U)  summarizes  the 
statistics  of  S in  the  limit  a = 0.  Special  attention  is  given  to 
statistics  in  frequency  and  pulse  propagation.  The  corrections  to  the  limiting 
statistics  are  derived  in  App.  (B)  and  discussed  in  Sec.  (5).  The  appearance 
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of  multiple  Fermat  paths  is  demonstrated  in  Sec.  (6).  Media  with  mviltlple 
scales  are  introduced  in  Sec.  (7)  &nd  the  distinction  between  full  and 
partial  saturation  is  made.  In  the  fully  saturated  case  there  is  a simple 
modification  of  the  results  for  a single  scale  medium  (Table  2).  The 
partially  saturated  regime  is  more  difficult.  Sec.  (8)  is  devoted  to  partial 
saturation  in  a medium  like  that  encountered  in  atmospheric  optics.  App.  (C) 
contains  some  calculations  relevant  to  Sec.  (8)  and  App.  (d)  discusses  some 
other  kinds  of  multiple  scale  media.  Methods  for  handling  inhomogeneous  and 
anisotropic  media  are  given  in  Sec.  (9)  and  Apps.  (E)  and  (F). 

2.  FIRST  .\ND  SECOND  MOMENTS  FROM  THE  PATH  INTEGRAL 

Feynman^  pointed  out  that  the  solution  to  Eq.  (l.U)  with  the  boundary 

condition  in  Eq.  (l.5)  is  given  by  an  infinite  dimensioneil  integral.  It  is 

defined  as  the  limit  of  a finite  dimensional  integral  with  2n  - 2 integration 

variables  corresponding  to  the  Cartesian  components  of  n - 1 two-dimensional 

vectors  r. , j = 1,  2,  . . . n - 1.  With  the  convention  that  r . I . „ and  r, I , 

J JJ=OJ'J=n 

are  the  source  r^  and  receiver  r coordinates,  Feynman's  integral  is 


C(r,T^,t)  = 


n 

r “ / 

/ kn  \ 

ikR  r /] 

^ L 1 jy 

[ 2iTiR  ) 

n , tV 

j\j=i  J/ 

(2.1) 


where  each  component  of  rj , J =1,  ...  n - 1,  is  integrated  over  the  range 

-00  to  +«  and  Zj  = JR/n.  In  p,  Tj  is  understood  to  be  a vector  in  the  (x,y) 

plane  and  e is  a unit  vector  in  the  z direction.  At  each  point  r,  , r„...r  , 
z 1*2  n-1 

in  the  integration  volume,  the  n-1  points  in  space  (r^,z^),  (r2,Z2),  ... 
(r^_^,z^_^)  can  be  thought  of  as  discrete  points  along  a path  r(z)  connecting 
'to  (r,R)  with  rj  = r(zj),  see  Fig.  (U).  In  this  sense  Feynman's 
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I 


r 


integral  is  an  integral  over  paths.  Associating  R/n  with  a differential 
increment  dz  in  range  the  argument  of  the  exponential  has  a continuum  limit 


dz 


<;^exp  l^-ik  I J u(?„(z)  + dz]  ^ 


R R 


= exp 


- V ^m^n j jp(/(^n^"^  " ’^-^m) 

0 0 


(2.M 
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(2.2) 

where  r'  = dr/dz.  The  path  integral  for  S can  then  he  schematically  written 
as 

r R 2 “I 

<?(r,r^,t)  = ^ / d(paths)  exp  ik  j -u^r(z)  + e^z.tjj  dz J 

(2.3) 

where  the  integration  is  over  all  paths  connecting  (r^,0)  to  (r,R)  and  the 
volume  element  in  path  space  d(paths)is  the  coefficient  of  the  exponential 
in  Eg.  (2.1) . 

We  will  be  computing  averages  of  products  of  path  integrals  and  the 
following  formiila  will  be  needed.  Let  r^{z),  n = 1,2,...  be  some  set  of  paths 

and  + 1 corresponding  phases.  Then  if  either  u is  a Gaussian  random  field 

✓ 2 1—6 
or  kL  / <<  1 and  its  statistics  are  arbitrary,  it  is  well  known  that 


I 


As  a first  application  of  the  path  integral  we  can  compute  {S) . 
This  will  not  turn  out  to  be  a particularly  interesting  quantity  but 
the  calculation  is  simple  and  it  will  show  how  path  integrals  work  and 
where  the  Markov  approximation  comes  in.  Bringing  the  average  inside 
the  path  integral  and  using  Eq.  (2.4)  yields 

RH 

rs 


<« . i 


d(paths)  exp 


+ (z  - 
(2.5) 


^,o) 


dzdz' 


The  Markov  approximation  now  appears  as  follows.  The  parabolic  wave 
equation  assumes  that  the  normals  to  the  wave  fronts  point  in  directions 
that  are  close  to  the  z-axis.  In  terms  of  the  path  integral  this  means 
that  for  the  important  paths  jr^l  = jdr/dz]  must  be  small.  It  then  follows 
that  for  important  paths  [ r(z)  - r(z' 

Eq.  (2.5)  becomes 


^r(z)  - r(z')j 


^ + (z  - z')^  « Iz  - z^l^and 


r 2 r /•”  1 

r R 

I-  — / #p(|z  - z l,0)dzdz'l  fd(paths)  exp 

— /(r  ) dz 

L 0 0 -1*' 

^ •/ 

^ 0 ^ 

(2.6) 


The  remaining  path  integral  is  Just  the  path  integral  for  5 and  for 

o 

R » L the  double  Integral  over  0 can  be  replaced  by  Rp(0,0).  The  final 
result  is  then 


<(J>  = <S  exp  [-J4  ] 
o 


(2.7) 


3-5 


This  is  the  usual  formula  obtained  in  the  Markov  approximation.''  ' What 
we  have  seen  here  is  that  this  approximation  has  a very  natural  inter- 
pretation in  terms  of  the  path  integral  and  that  it  is  valid  as  long  as 
the  parabolic  wave  equation  is  valid. 
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h 


If 


I 


Since  ^ is  large  in  the  saturated  region  (S)  is  exponentially 

small  and  therefore  not  particularly  interesting.  The  same  is  true  for 

<(S<5>  and  its  complex  conjugate  <4*5*).  The  path  Integral  for  <44>  will 

be  a double  path  integral  over  two  paths  r (z)  and  r (z)  and  will  con- 

1 2 

tain  a factor 


R R 


1 

(V 

%(.Z)  - r^(z')] 

+ (z  - z')  ,oj  + 2pl^^r^(z)  - + (z'  - z)^,t||dzdzj 

(2.8) 

where  t is  the  time  difference  between  the  two  <S ' s in  the  average.  This 
factor  is  of  order  exp  ] in  all  important  regions  of  path  space  and 
(i?S)  is  exponentially  small. 

A more  interesting  quantity  is  <4*(2)4(1)>  where  4(1)  is  a shorthand 
notation  for  4*(2)  for  formula  for  this 

object  is 


(2.9) 


where  the  path  integral  is  a "double  path  integral"  over  two  paths  r (z) 
and  r^(z)  connecting  (r^^j^.O)  to  (r^,R)  and  (r^^.O)  to  (r^.R)  respectively 
and 


R R 


= \ff\ ( + (z  - z')^, o) 

+ + (z  - z')^,o) 

- 2p(V^r^(^)  - r2(z')f  + (z  - z')^t^  - t^jjdzdz' 


(2.10) 
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> 


There  is  now  a region  in  path  space  where  the  integrand  is  not  exponen- 
tially  small.  It  is  r (z)  w r (z)  and  almost  all  of  the  path  integral 
will  come  from  this  region.  As  before,  (r^^Cz)  - r^(z  ) ) and 


/“>  / \2  / 2 
^r^(z)  - r^Cz  ) I can  be  neglected  relative  to  (z  - z ) and  in  the  same 


2 2 
spirit 


. 2 2 

(z)  - r (z')]  + (z'  - z)*^  can  be  approximated  by 


f — > 2 2 — 

a!  ( r (z)  - r + (z^  - z)  where  z = J(z  + z').  Then  for  R » L 

1 ' 1 2 ' 


the  integral  over  z - z^  can  be  done  and 

R 

/ 

V 

0 


(z)  - r^(z) I , t 


1 - S) 


dz 


(2.11) 


where 


= k^{p(0,0)  -p(lr|,t)j 


(2.12) 


At  this  pniJit  it  is  convenient  to  cliange  variables  to  paths  u(z)  and 


v(z) 


u(z)  = i^r^(z)  + r^(z)j  - 


v(z)  = r^c.".)  - r^Cz) 


^2^R  " + ^02> 


t-l) 


(2.13) 


wiiieh  satisfy  the  endpoint  conditions  u(0)  = u(r)  = ^ and  v(0)  = ’ 

v(R;  = (r^  - r^j . From  the  finite  form  of  the  path  integral  in  Eq.  (2.1)  it 
Is  clear  that  this  change  of  variables  is  allowed  and  that  the  associated 
Jacobian  is  equal  to  unity.  After  integrating  the  first  term  In  the 
exponential  by  parts  and  using  the  endpoint  conditions,  the  path  integral 
for  (2)4"(1)^  becomes 

R R 

(2)S(^)y‘  = ^d^(paths)  exp  |^-  Ik  ^u(z) -v" ( z)dz  - ^ d(|v(z)| 

(2.11*) 
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,t)d 


I 


where  t = In  analogy  with  the  formula 


hP 


dy  = 6(x) 


(2.15) 


the  integral  over  the  path  u in  Eq.  (2.14)  will  produce  a "6 -functional" 
which  forces  to  vanish  identically.^^  With  the  endpoint  conditions 


given  above  v(z)  must  then  be 


?(z)  . 


(‘-i) 


(2.16) 


In  d the  path  v can  then  be  replaced  by  the  right-hand  side  of  Eq.  (2.16) 
and  the  factor  containing  d then  becomes  just  exp  [-Jd].  The  remaining 
path  integral  is  equal  to  ( 2TiR/k )~^and  the  result^"^  reproduces  Eq.  (l.lU) 


<(S*(2)($(1)) 

5*(2)(S^(1) 


= exp  [-JD(1,2)] 


(1.14') 


where  D(l,2)  is  a shorthand  notation  for  D(r  - r„,  r - r , t - t„) . 

1 2 ol  o2  1 2 


Appendix  A contains  an  explicit  calculation  of  the  first  correction 

to  the  Markov  approximation  for  <(S*(2)(S (l)> . It  is  shown  to  be  propor- 

,2.  1 

tional  to  the  rms  multiple  scattering  angle  (\p  )R/L)  which  must  be 
small  if  the  parabolic  wave  equation  is  valid.  Henceforth,  all  calcula- 
tions will  be  done  in  this  Markov  approximation.  The  general  prescription 


is  that  whenever  p 


is  to  be  re 


never  p (^(tiCz)  - r^(z'))^  + (z  - z')^,t^  - tjj  appears,  it 
placed  by  6(z  - ==')P 


Turning  now  to  the  calculation  of  ((S*(u')(S(a))) , the  path  Integral 
for  this  quantity  will  contain  (with  k = k(a))  atid  k'  = k((jj')} 


n n 

^xp  [-ikypJr^Cz)  + e^z)dz  + ^^(z) 

N I 0 ***  0 


+ e z)dz  \ (2. : 

J/ 


Where  the  time  dependence  of  p has  been  suppressed  and  the  subscript 
indicates  that  for  a dispersive  medium  p can  depend  on  u.  Let  us  first 
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compute  this  average  in  the  absence  of  dispersion.  When  is  indepen- 
dent of  u it  is,  in  the  Markov  approximation 

R 

exp  -J(k  - k')\p(0,0)  -kk'/[p(0,0)  - p(lr  (z)  - r (z)|,0)]dz 

0 1 2 J 


(2.18) 


For  paths  which  make  a significant  contribution  to  the  path  integral, 


the  second  term  in  the  argument  of  the  exponential  must  be  of  order  unity 

/ -2 

or  less.  In  this  term  one  can  therefore  approximate  kk  by  k where 
k = i(k  + k').  Generalizing  to  dispersive  media,  one  finds  that  in  the 
same  approximation  the  result  is  just  Eq.  (2.18)  with  Rp(0,0)  replaced 
by  where  was  defined  in  Eq.  (l.20).  The  path  integral  will 

also  contain  a factor 


^r^(z)j  dz 


(2.19) 


which  can  be  simplified  by  making  an  orthagonal  transformation  to  paths 
u and  V defined  by 


k'v(z) 


->  / \ \ k'vU 

ri(z)  = u(z)  - 


-»■  / \ \ kv(z) 

r^(z)  = u(z)  - -1^ 


(2.20) 


After  making  this  transformation  the  path  integral  factors  into  a product  of 
integrals  over  u and  v.  Upon  dividing  by  the  integral,  over 

u cancels  and  the  final  result  is,  for  |u>  - a)'|  small  compared  to  o)  = ^ ((iH--m'), 


(o)’ ) ^ (m)y  _ 

(o)' )(?  (o)) 
o o 


i f 

: ^ V “g  ; j 


A ( 0)  - »i) ' ) 


(2.21) 


L 


where 


1 d( paths) 

ik  ^ 

CVJ 

o 

rr2 

1 [|5(0,0)-(S(v(z) 

2(k-k') 

- k 

0 

J 

► 

d( paths)  exp 

ik^  f 

2(k-k’)  J 
” /\ 

R 

(v'( 

z))  d^ 

0 


(2.22) 


In  the  saturated  region  where  $ is  large  | v(z) j will  be  very  small  for 
the  important  paths  and  the  expansion 


p(0,0) 


- p(l  v(x)  1 ,0)  as  ip(0,0) 


(2.23) 


can  be  used.  The  path  integral  for  A is  then 


(2.2h) 


This  type  of  path  integral  was  evaluated  by  Feynman  and  setting  k - k'  = 

(a)-a)')/c  it  is  equal  to 
6 


A((i)  - 0)') 


Li 

\ / 

/ (D  - U)'\i 

n(6i  — ) 


(2.25) 


-2  ^ 2 / h 3 

where  ojq  = Rp(0,0)  /c^  and  a = 6^L  /R  ^(0,0M  = H/i'.  Combining  equations 
(2.21)  and  (2.25)  yields  Eq.  (1.19).  Some  features  of  these  correlations 
in  w were  mentioned  in  the  Introduction.  We  will  return  to  their  inter- 
pretation in  Section  VI. 
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Except  for  the  explicit  verification  of  the  validity  of  the  Markov 
approximation,  the  above  results  could  be  obtained  by  more  familiar  techniques 
which  do  not  employ  the  path  integral.  The  power  of  the  path  integral  will 
become  apparent  in  the  next  section  when  higher  order  moments  are  computed. 
They  are  extremely  difficvilt  to  treat  by  the  usual  techniques. 


3.  HIGHER  MOMENTS  FOR  SMALL  a 


When  i is  large,  the  average  of  any  path  integral  will  be  exponen- 
tially small  unless  there  is  a region  of  path  space  where  each  path 

associated  with  an  (?  is  close  to  a path  associated  with  an  S*.  Such  a 

2 

region  does  not  exist  for  (5)  or  (S  ) and  we  have  already  seen  that  they 
are  exponentially  small.  More  generally,  any  moment  with  an  unequal 
number  of  (S’s  and  <5  's  will  be  vanishingly  small. 

Beyond  (5*(2)<S(1))  the  first  nontrivial  object  is  (5*(4)(S(3)(S*(2)(S(1)) 
It  is  given  by  the  quadruple  path  integral  over  four  paths  r^ (z) . . . r^ (z) 


((S*(4)(S(3)(S*(2)(S(1))  ={2k)~^  Jd  (paths)exp|-^ 


J=1  0 


(z))  ^dz 


(3.1) 


where 


4 . 

= ^ - rj(z)l,t.  - tj)dz  (3.2) 


i,J=l 


There  are  two  regions  of  path  space  where  M is  of  order  unity  or  smaller. 
They  are:  (a)  Ir  (z)  - r (z)]  < L/$ , | r (z)  - r (z)]  < L/$  with  the  dis- 
tance  between  pairs  of  paths  arbitrary  and  (b)  ir^(z)  - r^(z)l  < L/^ , 

I ^ I 

|r  (z)  - r (z)(  < L/$  again  with  the  distance  between  pairs  of  paths 
arbitrary.  In  region  (a)  where  Ir  (z)  - r (z) | is  of  order  L/^ , the 
oscillating  factor 


I 


exp 


dz 


exp 


(z)  + r''(z) 


(3.3) 


in  the  path  integral  will  restrict  Ir^^(z)  + r^^{z)l  to  be  of  order 

X M 

I ^ I 

2$/(kLR).  For  a typical  path  | r (z)  + r (z)|  will  then  be  roughly 

X 

— I — ) (r^^(z)  + r^^Cz))  ~ . The  centroid  of  the  other  pair 

will  be  restrained  in  a similar  way.  It  follows  that  most  paths  will  be 

such  that  the  ratio  of  the  distance  between  the  pairs  to  the  scale  length 

2 -1 

L is  roughly  #R/(6kL  ) = or  , where  a is  the  parameter  defined  in  the 
Introduction.  For  small  a the  pairs  are  separated  by  many  times  L and 
therefore  are  uncorrclated.  In  region  (a)  M then  reduces  to 

R / 

M /d(lr^(z)  - r^Cz^l.t^  - t2)dz  + - r^(z)l,t^  - t^)dz 

0 0 

, (3.4) 


and  in  region  (b)  it  becomes 

R 
i 

M 


n 


R 


^(z)  “ ^2^^^  * /*d(lr^(z)  - r^(z)|,t^ 

0 0 


t )dz 
4 

(3.5) 


Thus  in  each  of  the  two  important  regions  of  path  space,  the  quadruple 
path  integral  factors  into  the  product  of  two  double  path  integrals,  each 
of  which  is  precisely  the  integral  encountered  in  the  calculation  of 
((S*(S) . The  result  is  that 

<(S*(4)5(3)(S*(2)(S(1)>  « <<S*(4)(5(3)><(5*(2)(5(1)>  + <5*(2)<5(3)><(J*(4)5(1)> 

(3.6) 

where  the  two  terms  come  from  the  two  regions  (a)  and  (b) . 
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In  Appendix  B the  error  in  Eq.  (2.6)  is  obtained  by  computing  the  first 
correction.  It  is  of  order  a and  will  be  discussed  in  detail  in  Section  V. 


Generalizing  to  an  arbitrary  moment  is  easy.  The  general  non- 

vanishing  moment  Is^FI ‘S*(j)  ntJd^and  can  be  written  as  an  integral 
%1  i=l 

over  2n  paths  r (z)  and  r (z).  Tliere  will  now  be  nl  important  regions 
of  path  space  corresponding  to  the  number  of  ways  paths  r (z)  can  be 
paired  with  the  paths  r.(z).  In  each  of  these  regions  the  2n-tuple 
path  integral  can  be  approximated  by  a product  of  n double  path  integrals. 
Some  simple  combinatorics  shows  that  the  result  will  be  as  follows.  Let 
i be  a permutation  of  the  indices  i.  For  example,  if  n = 3 and  the 
permutation  is  (1,2,3)  -»  (3,1,2)  then  1=3,  2=1  and  3 = 2 or  if  the 
permutation  is  (1,2,3)  ->  (2,1,3)  then  1=2,  2=1  and  3=3.  With  this 


notation 


/ n <s*(3)  n 5(^«  n <^*(3)c?(i^ 

^j=l  i=l  ' perms  i,j=l 


(3.7) 


wliere  the  sum  of  over  all  n.'  possible  permutations  of  the  indices  i. 


Tlie  same  result  holds  for  correlations  in  frequency.  Extending  the 

notation  iS(j)  to  include  a frequency  label  m we  have 

j 


= exp  -iD(i,J)  - m ) 

o o \ g / J 


(3.8) 


which  liolds  in  the  saturated  region  where  only  small  values  of  u ~ w , 

|r^  - ^ j 1 I etc.  are  interesting.  The  same  construction  that  led  to 
Eq.  (3.7)  for  equal  frequencies  then  shows  that  it  holds  for  unequal 
frequencies  as  well. 

The  interpretation  of  Eq.  (3.7)  will  be  given  in  the  next  section. 

A final  remark  here  is  that  the  arguments  loading  to  Eq.  (3.7)  do  not 
depend  on  the  validity  of  the  Markov  approximation.  The  latter  is  needed 
only  when  (l)^(2)^  is  explicitly  evaluated. 


U.  tie  statistics  of  (S  in  the  limit  a = 0 


The  moments  of  Eq.  (3.T)  correspond  to  a complex  Gaussian 

distribution.  The  probability  that  5(j)/5  (j)  will  he  equal  to  11.  for 

o J 

j = 1, . , . n is  then 


n 

^n^^l’'"^n^  = |det[TrM  ]j  ^exp  ^i^'^^ij 


(4.1) 


where  the  n by  n matrix  M is 


(0)  - u \2 


Eq.  (1.10)  corresponds  to  the  special  case  n = 1 and  Eq.  (1.18)  is 
obtained  by  dividing  P.jCil.Tl')  by  P^(Tl).  The  measure  is  d^q  = d(lmn)d(Ren) . 

In  principle,  Eqs.  ()t.l)  and  (1*.2)  determine  all  the  statistical 
properties  of  d . For  example,  it  follows  from  Gaussian  statistics  that  for 

<S(j)  = A(j)e^‘*’  the  correlations  of  amplitude  and  rate  of  phase 

A - di_  18 

(ft  = — ^ 


~^(l<?y  )>  = E(!m^2|)  - hil- 

where  E and  K are  the  complete  elliptic  integrals  of  the  first  and  second 


kinds  and 


<(P(1)9(2)>  = - £n(l-  |M^g|  ) 


(4.4) 


where  in  the  secono  line  the  expansion  of  D in  Eq.  (1.17)  has  been  used. 

Eq.  (4.4)  can  bo  extended  to  the  derivatives  of  cp  with  respect  to  r and 

-»  , .2 

r^  in  the  obvious  way.  Intensity  correlations  l(j)  = |A(j) | are 


22 


simpler  with 


<I(1)I(2)>  = (I)^  + 1<(S*(1)5(2)>1^ 


(4.5) 


The  appearance  of  Gaussian  statistics  in  frequency  is  somewhat 
unf ariiiliaj’.  To  see  what  it  implies,  let  us  compare  the  saturated  regime 
to  a simplier  unsaturated  one.  When  simple  first  order  geometric  optics 

„ ilDT 

applies  and  the  medium  is  nondispersive,  5(uj)  is  equal  to  o (uj)e  where 

o 

T is  a fluctuating  time  shift  independent  of  m.  Under  such  propagation 
conditions,  the  statistics  in  cu  are  essentially  trivial.  The  envelope 
cS(uo)  fluctuates  but  does  so  in  such  a v;ay  that  at  a fixed  time  when  t 
has  a definite  value  a knowledge  of  S at  one  value  of  aj  determines  S for 
all  u).  Another  way  to  say  the  same  thing  is  that  a pulse  will  be  sub- 
jected to  a random  time  shift  but  will  not  be  distorted  in  shape.  For 
propagation  in  the  saturated  regime  tlie  statistics  of  5(uj)  are  non- 
tj'lvial  and  things  are  completely  different.  At  one  fixed  time  a knowledge 
of  <S (u))  at  one  u)  yields  only  statistical  infoj’mation  about  3 at  nearby  fre- 
quencies. Correspondingly,  the  medium  will  distort  a pulse  in  a way  that 
is  predictable  only  statistically.  A peculiarity  is  that  (5* (o) ')d  (ai)) 
and  has  a phase  corresponding  to  an  average  retardation. 

Tlie  above  remarks  about  5(a))  are  most  easily  made  quantitative  in 
terms  of  pulse  propagation.  It  is  worth  going  into  this  in  some  detail 
both  because  tlie  physics  is  interesting  and  because  it  will  connect  with 
the  Fermat  paths  of  Sec.  (6).  For  simplicity  the  unperturbed 
medium  will  be  assumed  to  be  nondispersive  with  w = ck.  Let  the 

e f ('ju)duj  whore  f (-uj)  = f (in), 
o o o 

Tal-.ing  the  unperturbed  arrival  time  as  tlie  origin,  the  received  signal 

/-iujT  . ~ 

e 0(uj)f  (cu)du).  The  signal  f (t)  is  a Gaussian 
o r 

random  variable  whoso  complete  statistics  are  determined  by  the  covariance 

~ 19  ~ 

of  f 5.  Assuming  that  variations  in  S (iu)f  (oj)  over  a frequency 
o o o 
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corresponding  to  the  width  of  (<S*(-iu))5 ( can  be  neglected,  this  co- 


! 


variance  is 


<f*(u)')<S*(x')f  (u))(S(u))>  = 1?  (J(“>  + “>'))  1*  - “>') 


(J*.6) 


with 


where  the  small  a limit  of  Eq.  (1.19)  has  been  used  and  = a)^a/6  = 

CgL  /R^p(0,0)  . 

2 

Denoting  the  received  intensity  f (t ) by  the  average 

r 

is  a measure  of  the  distribution  of  energy  over  arrival  times.  Accord- 
ing to  Eq.  (4.6)  it  is 


^ (t  )<iT  y 


00 

*/ 


e A(u))da)  = 


«.  - 2 2 
« r / \n/  \2  -n  n u:,  t 

-2w^  I (-1)  (ntr)  e 1 

- n=l 


T < 0 


T > 0 


(»*.7) 


and  vanishes  for  t < 0 because  the  integrand  is  analytic  in  the  upper 
half  plane.  Evidently,  all  the  energy  comes  in  after  the  unperturbed 

< < r 2^-1 

arrival  time  and  is  confined  to  a region  0 «.  t ~ ('u^n  ) . Hie  net 

retardation  is  consistent  with  w>.at  was  said  above  about  the  phase  of 
{S*  (m)) . The  complete  abse  ce  of  energy  for  t < 0 is  peculiar  to 

the  limit  of  small  a and  will  later  be  seen  to  have  a simple  physical 
interpretation  (see  Sec.  (6)). 

For  a sharp  transmitted  pulse  the  distribution  of  energy  over  arrival 
times  can  be  thought  of  as  being  due  to  two  effects.  One  is  the  wander  in 
arrival  time  of  the  center  of  the  pulse  and  the  other  is  spreading  of  the 
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pulse  around  Its  center.  The  two  effects  are  In  principle  distinct, 
the  simple  case  of  propagation  in  an  unsaturated  regime  where  (S(uj)  = e 
the  wander  is  of  order  (t  ) while  the  spread  is  just  the  width  of  the 
transmitted  pulse.  As  we  will  now  see,  in  the  saturated  regime  the 
spread  and  wander  are  roughly  equal.  The  width  of  measures  the 

sum  of  spread  and  wander.  A quantity  which  measures  the  spreading, 
independent  of  wander,  is 


For 

1(UT 


<y^(T  + T 


")^(T')dT  '> 


P(t)  = 


00 

^ J ^J(T'))dT'y 


(4.8) 


When  f 5 has  a Gaussian  distribution  P(t)  is 
o 


P(T)  = 


P (T) 
O 

2-!i 


CD  CD 

/lA(ca)l^da)  + ^fe 


iiuT . ,2 

lA(()i))  da) 


(4.9) 


where 


P (t)  = 
o 


i f V 

\ y f (T  + T ')f  (T  ')dT  'I 

'-CO  O Of 

(i 


f^(T')dT'r 


(4.10) 


Tlie  two  terms  in  P(t)  have  the  same  height  at  t =0.  The  spike  propor- 
tional to  P (t)  falls  rapidly  leaving  the  second  term  whose  width  is  a 
o 

measure  of  the  spread.  Comparing  with  Eq.  (4.7)  one  sees  that 


00 


1a(u))1 


2 

du) 


00 


'<JJ(T  + T')>(JI,(T')>dT' 


-CO 


(4.11) 


and  it  is  clear  that  the  spread  and  wander  are  essentially  the  same.  A 
physical  Interpretation  of  the  two  pieces  of  P(t)  will  be  given  in 
Sec-  (5).  Finally,  a useful  formula  is 
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(4.12) 


|A(id)1  = 


sin 


w 


2/  u)  U 
+ sinh  /- — r 

\2u>i; 


It  is  interesting  to  ask  why  it  is  that  the  square  of  the  auto- 

2 

correlation  of  f rather  than  the  autocorrelation  of  f appears  in  P . 

o o o 

The  answer  is  that  when  f (u))<S(u))  has  a Gaussian  distribution,  the  medium 

o 

cannot  transmit  any  information  that  is  not  contained  in  the  coherence 

(f*(iu^)(S*(u)^)f  (a))(S(u))).  As  given  by  Eq,  (U.6  ) this  coherence  depends 
o ^ o 

only  on  jfCu))!^  and  the  medium  can  only  transmit  information  about  the 
autocorrelation  of  f. 

The  statistics  of  the  signal  as  a function  of  spaclal  wave  numbers 

can  be  analyzed  in  a similar  way.  Multiplying  o(r,r  , t)  by  a suitable 

o 

function  of  r and  integrating  over  r one  can  represent  a boundary  con- 
o o 

dition  at  z = 0 corresponding  to, say,  a plane  wave  emerging  from  a finite 
aperture.  For  such  a signal,  the  Fourier  transform 


5 (f)  = 

o 


(2it) 


-ir  • 2 


S (r)d  r 


(4.13) 


will  be  sharply  peaked  around  some  Z = I , With  the  correspondences 

o 

t -*  l,  u)  r one  can  proceed  as  above  and  discuss  spread  and  wander  in 


f.  Again  the  medium  can  only  transmit  information  contained  in 
(jS* (,r)S  (r  *))  which  will  typically  depend  only  on  |(S  (J(r  + r^))|  . 


5.  CORRECTION  TO  THE  a = 0 LIMIT 

The  leading  corrections  to  the  0=0  limit  are  computed  in  Appendix  B. 
The  main  results  are  as  follows. 

The  order  a correction  to  (l**)  is  dominated  by  fluctuations  near  the 
transmitter  and  receiver.  This  is  not  unexpected  since  near  their  end- 
points the  paths  cannot  be  separated  into  uncorrelated  pairs.  Explicitly 
<l">  is 


26 


(l°)  = + Jn(n  - l)aC  + 0(ck^)) 


(5.1) 


where 


_ (3ti)^ 


/'2~, 

q p(q,0)dq 

OO 

/qp  (q,0)dq 


(5.2) 


and  p is  the  three-dimensional  Fourier  transform  of  p 


w 

i(|xl,t)  = stn(q|x|)p(q,t)dq 


(5,3) 


The  consequences  of  the  fact  that  the  error  grows  with  n were  noted  in 
the  Introduction.  Note  that  the  correction  to  (l*')  is  positive.  This 
means  that  the  intensity  fluctuations  overshoot  (i.e.,  become  larger 
than  Rayleigh)  near  the  boundaries  of  the  saturated  regime. 

The  correction  to  a general  correlation  can  also  be  computed.  They 
are  always  fractionally  small.  For  example,  (5* (2)5* (2)4 (1)5 (1))  is 

. , , “D(l,2)  rt  . • • ..u  — .4- 


proportional  e 

^ -D(l,2) 

order  ae 


in  the  a = 0 limit  and  the  correction  to  it  is  of 
, The  corrections  to  Intensity  correlations  are  the  most 


interesting.  In  the  a = 0 limit  (l(t  )I(t  ))  is  equal  to 

X M 

, , 2 f -D(t  - t ) 

(I)  1+e  1 2.Att=t  the  order  a correction  is  given  by 

I J 1 2 2 

Eq.  (5.1).  However,  at  large  |t  - t [,  (l(t  )I(t  ))  must  approach  <I) 

1 ^ 1 45 

and  the  corrections  must  go  to  zero.  It  turns  out  that  half  the  correc- 
-D(t  - t ) 


tion  dies  like  e 


2 but  the  other  half  falls  much  more  slowly. 


leading  to  a coherence  tail.  (Note  that  this  is  consistent  with  what 
was  said  above  about  the  corrections  always  being  fractionally  small.) 
For  the  general  intensity  correlation  the  coherence  tail  is 


(5.4) 


, ,2  -D(l 

<1>  1+e 


r 

,2)  «/r,  % 


fqP(q,0)dq 
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6.  FERMAT  PATHS 

There  is  an  interesting  connection  between  averages  of  the  path 
integral  and  averages  over  Fermat  paths  which  satisfy  the  perturbed 
ray  equation 


r'  (z)  + 7p(r(z)  + e^z)  =0 


(6.1) 

fh  b\ 

where  V = — ir~  • This  will  be  Illustrated  for  the  special  case  of 
Vax  hyj 

soiirces  and  receivers  located  at  r = r = 0 so  that  4 is  a function  only 

o 

of  time.  The  path  integral  for  (4*(t)4(0))  is 

R 


<4*(t)4(0)>^^^  yd^(paths) 


exp 


iky*[i(r'(z))^  - - p(r^(z)  + e^z,o) 


+ p\r  (z)  + e z 


,t) 


dz 


(6.2) 


In  the  saturated  region  we  know  that  for  ((5*(t)(S(0))  to  be  non- 

I “>  ->  I 

vanishing,  t must  be  small  and  that  only  paths  for  which  |r  (z)  - r (z) | 

^ 1 ^ 2 ^ . 

is  small  (~  h/i)  contribute.  Changing  variables  to  w(z)  = jlr  (z)  + r (z)) 

'X  2 

and  v(z)  = r (z)  - r (z),  the  path  integral  can  then  be  approximated  by 
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‘(t)(5(0))=^^yd^(paths)  exp  -Ikyj  v(z)  w"(z)  + Vp  w(z)  +e^z,o|| 


- tp(w(z)  + e^z,o)|dz  > (6.3) 


where  the  first  term  In  the  argument  of  the  exponential  has  been  integrated 

by  parts  and  a dot  indicates  differentiation  with  respect  to  time.  The 

integration  over  the  path  v(z)  produces  a "6 -functional"  which  forces 

w + Vp  to  vanish  for  all  z.  Thus  the  Integral  over  w(z)  is  restricted 

to  paths  which  satisfy  the  ray  equation  (6.1).  This  is  a general  feature 

of  the  saturated  region.  Higher  order  correlations  are  dominated  by  con- 

figurations  where  paths  r.(z)  and  r ,(z)  are  pair-wise  close.  A similar 

analysis  shows  that  for  each  such  pair,  the  path  w (z)  = i r (z)  + r ,(.z) 

i I i 1 J 

satisfies  Eq.  (6.1). 

Equation  (6,3)  can  be  further  analyzed.  For  most  media  p and  p are 

statistically  independent.  The  average,  (),  can  then  be  thought  of  as 

two  independent  averages  ( ) and  ( )•  over  p and  p.  The  6 -function  of 
^ p p ^ 

v'  + Vp  that  is  produced  by  the  Integration  over  v is  effected  only  by 
the  average  ( ) while  the  phase  exp[ikt/p]  is  effected  only  by  the  other 
average  ( )*.  It  is  therefore  possible  to  write  (<S*(t)5(0))  as  the 
integral  over  paths  w of  a p-averaged  6 -functional  which  can  be  inter- 
preted  as  the  probability  that  a given  path  w will  satisfy  Eq.  (6.1) 
times  a phase  which  is  to  be  averaged  over  p.  To  do  this  correctly,  it 
is  necessary  to  go  back  to  the  definition  of  the  path  integral  in  Eq.  (2.1), 
The  integration  variables  vj^  and  wj^,  ^ = l,2,...,n  are  then  discreet  and 
the  mathematics  is  straightforward.  The  integration  over  the  v 's  can 
be  done  trivially  and  after  some  manipulation,  one  finds 

1^ 

(i5*(t)<S(0)>  = ^^d (paths) 'p(path)^xpj  ikt^p(w(z)  +e^z,o)dz|\  (6.4) 

. 

where  the  integration  d(paths)  is  over  paths  w(z>  With  a modified 
volume  element  * 
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d (paths) ^ = 


(6.5) 


which  does  not  contain  and 


, , /n\2n-2 

P(path)  = f- 


/n\2n-2  / 2/n  ^ 

(“)  \j".  ■ '•j’  * ’“‘"j  ‘ 


(6.6) 


is  the  probability  that  w will  satisfy  the  finite  difference  approximation 
2 


“r(w  ,+w.,-2w)+7n(w  +ez,0)=0 

jj2  j+1  j-1  j j z J 


(6.7) 


to  the  ray  equation.  In  the  limit  n -»  “>,P(path)  is  the  probability  (with 
a measure  d(paths)  ) that  w will  satisfy  Eq.  (6.1).  Equation  (6.4)  shows 
explicitly  that  (<S*(t)<S(0))  is  a sum  over  Eermat  paths  with  fluctuating 
phases  Finally,  bringing  the  average  of  p inside  the  exponential 

yields 


^<S*(t)S(0))  = ^^d (paths)  'P(path)sxp 


Vi 

(6.8) 

In  the  Markov  approximation  where  w(z)  is  neglected  relative  to  e 2 in 

2 

the  average  of  p,  Eq.  (6.8)  becomes 


<(S*(t)5(0)>  = exp 


] /d 


(paths) 'P(path) 


= exp  (-)  (I> 


(^rl< 


(6.9) 


which  is  the  standard  result. 


This  provides  a new  way  to  look  at  the  Markov  approximation.  It 
requires  that  an  average  like  <|^p(w(z)  + e z,o)dz|  > along  a path  which 
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satisfies  the  perturbed  ray  equation  (6.1)  should  be  well  approximated  by 

/ ^ \ 2 

the  corresponding  average  ((y'pCe  z,0)dz)  > along  the  unperturbed  ray. 

'0  ^ 

For  a homogeneous  and  isotropic  medium,  this  will  be  the  case  as  long 

2 i 

as  the  rms  multiple  scattering  angle  ((p  )R/L)^  is  small. 

According  to  Eq.  (6.8),  ^<S*(t)5(0))  can  in  principle  be  computed 

by  a geometric  optics  method  which  searches  out  the  rays  which  satisfy 

the  perturbed  ray  equation.  Geometric  optics  corresponds  to  an  approxi- 

6 

mate  evaluation  of  the  path  integral  by  the  method  of  stationary  phase. 

In  the  saturated  region  the  stationary  phase  approximation  will  in  fact 


be  valid  since  for  $ > 1 the  phase  ^ 


y^z  i 


s necessarily  large.  To  get 


Gaussian  statistics  for  5,  It  is  necessary  that  there  be  several  rays 

connecting  a given  source  and  receiver.  In  path  integral  language  this 

means  that  there  will  be  multiple  stationary  phase  points  and  S will  be 

a discreet  sum  EAue^'^^*  over  contributions,  one  from  each  stationary 
k 

phase  point  or  ray.  The  phases  and  amplitudes  A^^  as  well  as  the 
number  of  rays  will  fluctuate  with  p yielding  Gaussian  statistics  for  5. 

It  is  difficult  to  prove  rigorously  that  there  are  always  multiple 
rays  in  the  saturated  regime.  However  there  is  a simple  construction 
which  shows  the  essential  physics.  At  one  fixed  time  the  rays  are  sta- 
tionary points  of  the  path  length  S defined  by 


S = kf^'  [-i(r'(z))^  - u (r(z)  + e z)]  dz  . 
o z 


(6.10) 


Let  S(r)  be  S evaluated  for  the  special  paths  that  go  in  a straight  line 
from  the  source  at  (0,0)  to  an  arbitrary  point  (r,  z^)  with  0 < Zq<  R 
and  then  follow  another  straight  line  from  (r.z^)  to  the  receiver  at 
($,  R) . Multiple  stationary  points  of  S(r)  as  a function  of  r will  be 
indicative  of  multiple  stationary  points  of  the  complete  functional  in 
Eq.  (6.10)  and  the  spacing  of  such  points  in  r will  be  similar  to  the 
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1 


spacing  between  multiple  rays.  Now  doing  a simple  integral  shows  that 
S(r)  can  be  written  as 

S(^)  = h (?)^B  - S^(h  (6.11) 

kR 

where  B = — r- ^ and  S,  Is  kp  Integrated  along  the  above  mentioned 

path.  To  simplify  S(r),  B can  be  replaced  by  its  average  value  R/(6K). 
Then  defining  r = L u and  Sj^(r)  = $ f(u)  the  quantity  to  be  studied  is 

-*.2  -> 

fi  u - <t>£ (u)  (6,12) 

and  we  are  interested  in  its  stationary  points  which  satisfy 

Z -1>V  f(^)  = 0 . (6.13) 

By  construction  f is  a random  function  of  order  unity  which  changes  by 
order  one  when  its  argument  changes  by  order  one;  i.e.|Vf|~l  and  Vf 
changes  sign  roughly  each  unit  in  u.  For  >> 4>  the  first  term  in  Eq. 
(6.13)  dominates  and  there  will  be  a single  solution  near  u = 0.  This 
is  the  unsaturated  regime.  In  the  saturated  regime,  if  >>  f),  the  random 
I character  of  f guarantees  that  there  will  generally  be  many  solutions, 

spaced  by  about  one  unit  in  u (a  distance  L in  r)  and  filling  up  the 
Interval  0<|u|<  (0<  |r|  < ifL/n) . To  find  the  other  boundary  of  the 

saturated  regime,  41  > 1,  we  have  to  ask  when  the  multiple  rays  are 
physically  meaningful.  From  their  interpretation  as  stationary  phase 
points  of  the  path  Integral  it  can  be  verified  that  two  rays  will  be 
physically  distinct  if  S varies  by  a quarter  cycle,  i.e.,  order  unity 
between  the  two.  The  variation  in  S between  two  solutions  of  Eq.  (6.13) 
will  be  roughly  4>  and  if  they  are  to  represent  physically  distinct  rays 
4*  must  be  greater  than  unity. 

An  experiment  with  a pulsed  source  will  tend  to  see  several  arrivals 
corresponding  to  the  multiple  Fermat  paths.  This  random  multipathlng 
is  the  origin  of  the  rapid  fall-off  of  frequency  coherence  which  takes 
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place  In  the  saturated  regime.  To  see  how  the  orders  of  magnitude  work, 
the  difference  In  travel  time  between  the  ray  nearest  u = 3 and  the  fur- 
thest  one  out  at  |u  Is  (♦  # /(ZOm)  where  the  two 

terms  come  from  the  two  terms  In  Eq.  (6.12)  and  It  has  been  assumed  that 
$ >>  n.  Frequencies  which  differ  by  more  than  t^  ^ will  then  be  Incoherent, 
In  agreement  with  Eqs.  (1.19)  and  (1.21),  Note  that  t^  is  positive. 

This  Is  why  in  the  limit  a = 0 all  the  energy  arrives  after  the  unperturbed 
arrival  time  and  (t)  > vanishes  for  t < 0.  Also  the  two  terms  In  P(t) 
(Sec.  4)  can  easily  be  Interpreted  in  terms  of  fluctuating  multipath. 

The  spike  Pq(t)  is  the  autocorrelation  of  each  arrival  with  itself  and 
the  broad  second  term  is  the  autocorrelation  of  different  arrivals. 

Finally  a word  of  caution.  The  above  construction  vastly  underestimates 
the  number  of  rays.  In  reality  the  number  of  rays  is  probably  an 
exponential  of  <t>/fi  rather  than  di/n  as  the  construction  would  imply.  It 
may  be  extremely  difficult  to  actually  resolve  the  arrivals. 

It  is  interesting  to  consider  the  transition  into  the  saturated 
regime  in  terms  of  propagation  of  a pulse.  Consider  first  crossing  the 
line  $ = n from  the  region  where  both  $ and  12  are  large  but  12  < if. 

With  (1>  and  12  large  but  well  outside  the  saturated  region,  one  knows  from 
the  Rytov  approximation  that  the  receiver  will  see  a single  arrival  with 
a considerable  wander  in  time  of  arrival.  At  the  boundary  of  the 
saturated  region  the  pulse  will  begin  to  split  into  several  arrivals  and 
well  inside  the  saturated  region  there  will  be  many  arrivals  that  are 
spread  out  over  a time  long  compared  to  the  original  wander  in  the  single 
pulse.  Crossing  the  boundary  4>  = 1 from  the  region  where  both  and  12 
are  small  is  rather  different.  In  this  case  one  knows  that  well  outside 
the  saturated  region,  there  will  be  a single  arrival  with  no  discernible 
wander  in  time  of  arrival  accompanied  by  a small  scattered  wave  spread 
over  a continuum  of  arrival  times.  As  the  boundary  of  the  saturated 
region  is  approched , the  single  peak  will  shrink  and  the  scattered  wave 

t will  grow  in  amplitude.  Well  inside  the  saturated  region,  the  original 

peak  will  have  disappeared  completely  and  the  now  large  scattered  wave 
will  have  broken  up  into  a number  of  discrete  arrivals. 

• 
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7 . Media  with  Multiple  Scales 


So  far  it  has  been  assumed  that  the  fluctuations  in  u can  be 
characterized  by  a single  scale  size  L.  Technically,  this  requires 
that  the  expansion  of  0 , 

k^Ra  (1  r|  ,0)  = f 1 - ^ + a + ^ (7.1) 

\ 2L^  4L^  / 

-V  4 

through  order  r exists  and  that  the  coefficient  a is  of  order 
unity.  There  are  cases  of  practical  importance  where  this  is  not 
true.  For  example,  optical  index  of  refraction  fluctuations 
induced  by  Kolomogorov  turbulence  have  the  property  that  the 
(three  dimensional)  Fourier  transform  p (q)  of  p behaves  like 
|q(  over  a long  interval  in  q and  the  expansion  in  Eq.(7.1) 

makes  sense  only  when  the  cutoff  (inner  scale)  is  taken  into 
account  and  then  a is  very  large.  This  and  the  following  section 
are  devoted  to  these  media  with  multiple  scales.  It  will  be 
assumed  that  p (q)  goes  like  (q|  ^ ^ for  large  q where  4 > p > 1. 

(If  p is  greater  than  four  the  medium  acts  like  one  with  a single 
scale  size  and  for  p < 1 it  is  so  singular  that  does  not 

exist.)  In  practice  there  is  always  some  physical  cutoff  at  large  q 
(inner  scale) . However,  the  effects  of  such  a cutoff  will  be 
ignored  in  what  follows. 

For  p > 2,  the  length  parameter  L will  be  defined  by  Eq.(1.13) 
as  before  and  in  the  case  p < 2,  L will  be  defined  by 

M|r|  ,0)  = 0(0,0) 


(7.2) 
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for  small  Irj.  For  Kolomogorov  turbulence,  p is  equal  to  5/3 
and  p(0,0)  and  L are  related  to  Tatarskii's  by  2.91C^  = 0(0, 0)L 
The  parameters  fl*  and  n continue  to  be  defined  by  Eqs.  (1.6)  and  (1.7). 

The  main  qualitative  difference  between  propagation  in  single 
and  multiple  scale  media  is  that  in  the  latter  case  there  is  more 
than  one  saturated  regime.  In  terms  of  the  Fermat  paths  of  the 
last  section,  it  turns  out  that  in  a multiple  scale  medium  the 
smaller  scale  inhomogenieties  can  make  multiple  Fermat  paths 
before  the  large  ones  do.  This  leads  to  a new  kind  of  saturated 
regime.  Even  in  a single  scale  medium  with  p > 4 the  line  4>  = 

' 1 

is  not  a sharp  boundary.  In  reality  there  is  a transition  zone  | 

where  random  focusing  along  single  Fermat  paths  produces  intensity  i| 

M 

fluctuations  bigger  than  Rayleigh.  As  p decreases  below  four  this  s 

transition  zone  opens  up  and  becomes  a new  saturated  regime.  The  *| 

boundaries  of  this  new  regime  can  be  found  by  studying  the  object 
i fiu^  - $f (u)  of  Eq.  (6.12)  . 

To  see  when  the  smaller  scales  can  make  multiple  Fermat  paths,  > 


imagine  throwing  out  all  scale  sizes  larger  than  AL  where  1 > X > 0. 

The  new  scale  length  will  be  XL  and  $ and  f)  will  be  replaced  by 

X^‘^^4>  and  X^fi.  The  combination  i/Q  becomes  X and  is 

equal  to  unity  when  X = (<I>/n)^'^^^  . Thus  if  p < 4 the  small 

scales  can  make  multiple  Fermat  paths  when  4>  < fi,  i.e.,  before 
the  large  ones  do  at  4'  = n.  However,  if  these  multiple  paths  are 
to  be  physically  meaningful  must  be  greater  than  unity  and 

the  smallest  permissible  value  of  X is  <I>  Putting  everything 

together,  the  small  scales  can  make  meaningful  multiple  Fermat 
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paths  when  > 1.  This  is  one  boundary  of  the  new  saturated 

regime.  To  find  the  other  boundary,  we  need  to  ask  when  the 
multiple  Fermat  paths  can  be  separated  by  L.  For  a given  X the 
minima  of  ^X^S^u^  - X^^^$f(u)  extend  out  to  a maximum  |uj  which 
IS  the  largest  value  of  |u|  for  which  the  equation 
X^fiu  - (u)  can  be  solved.  The  maximum  |u|  is 

(iI>/£i)X^^  4)/2  noting  that  u is  distance  in  units  of  XL  one 

sees  that  the  Fermat  paths  can  be  separated  by  L when 
(<&/n)  X = 1.  For  p > 2 the  most  separated  paths  are  due  to 

large  scales  with  X = 1 and  the  other  boundary  of  the  new  region 
is  O = n.  However,  if  p < 2 the  smaller  scales  produce  the 
largest  separation  and  taking  the  smallest  permissible  value  <I> 
for  X one  sees  that  there  can  be  Fermat  paths  separated  by  L when 
> 1.  The  regime  where  there  are  meaningful  multiple 
Fermat  paths  all  lying  within  L of  each  other  will  be  called  the 
partially  saturated  regime . The  regime  where  the  spacing  between 
Fermat  paths  can  be  greater  than  L is  analagous  to  the  saturated 
regime  of  the  single  scale  case  and  will  be  called  the  fully 
saturated  regime  . The  boundaries  of  these  regimes  are  summarized 
in  Table  1. 


Partially  Saturated  Regime 

Fully  Saturated  Regime 

<I>>1,  4>/n>i 

$>i, 

Table  1.  Boundaries  of  the  Saturated  Regimes 


Although  these  boundaries  have  been  obtained  with  a hueristic 


Fermat  path  argument  they  are  in  agreement  with  what  one  finds 
from  more  precise  calculations.  It  is  known  that  outside  the 
saturated  regimes  the  intensity  fluctuations  ^ ) / <^ 

are  small,  implying  both  the  validity  of  the  Rytov  approximation 
and  the  absence  of  saturation.  Inside  the  saturated  regimes  (as 
given  by  Table  1)  the  intensity  fluctuations  as  computed  in  the 
Rytov  approximation  are  large,  signaling  the  onset  of  saturation. 
The  line  between  the  fully  and  partially  saturated  regimes 
corresponds  to  the  place  where  two  pairs  of  paths,  in  the  sense  of 
Sec.  (3) , can  be  separated  by  more  than  L.  When  they  are 
separated  by  more  than  L the  pairs  of  paths  are  completely 
independent  (full  saturation)  and  Gaussian  statistics  for  ^ follows 
immediately.  If  all  pairs  are  within  L of  each  other  (partial 
saturation)  then  one  expects  that  at  least  some  statistics  will 
not  be  Gaussian. 

Nothing  that  was  done  in  Sec.  (2)  or  App.  (A)  depended  in 
any  essential  way  on  the  assumption  of  a single  scale.  The  reader 
can  verify  that  Eg.  (1.14)  for  <^£*  (2)£’(1)  ^ at  equal  frequencies 
continues  to  hold  whenever  the  parabolic  wave  equation  is  valid. 

The  only  subtle  point  is  that  for  p < 2 the  rms  scattering  angle 
is  not  well  defined  and,  correspondingly,  in  App.  A,Eq.(A.6) 
cannot  be  approximated  by  Eq.(A.9).  However,  a rather  straight- 
forward analysis  of  Eq.(A.6)  shows  that  the  fractional  error  in 
Eq.(1.14)  is  of  order  D(k  ^,0)  and  it  is  known  that  D (k  ^,0)  < 1 
is  the  validity  condition  for  the  parabolic  wave  equation  when 
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p < 2.  Turning  to  coherences  in  frequency,  there  is  however  a 
significant  defect  in  the  theory  if  p < 2.  When  p is  less  than  2, 
the  path  integral  in  Eq. (2.22)  cannot  be  approximated  by  that  in 


.(2.24)  and  A must  be  understood  as  a function  defined  b 


Eq.(2.22)  whose  evaluation  would  require  a numerical  calculation. 


In  the  fully  saturated  regime  where  pairs  of  paths  can  be 
separated  by  L or  greater,  the  arguments  of  Sec.  (3)  proceed  as 
before.  One  readily  verifies  that  in  the  fully  saturated  regime 
the  statistics  of  £ are  Gaussian  and  the  discussion  of  Sec.  (4) 
applies  (except  Eq.(4.7)  which  assumes  Eq.(1.21)  for  A), 

Eqs.(B.12)  and  (B.17)  of  App.  B hold  in  the  multiple  scale  case. 

The  reader  can  then  verify  that  for  p > 2,  Eqs . (5 . 1) , (5 . 2) , (5. 4)  and 
(5.5)  for  the  corrections  to  Gaussian  statistics  continue  to  hold 
in  the  fully  saturated  regime  and  that  for  p < 2 these  same 
equations  hold  if  a is  replaced  by  a'  where 


- 4(p+l)^^Pr(3/p)  n 

" ^iS-2p)/p 


(7.3) 


The  situation  for  the  fully  saturated  regime  is  summarized  in 
Table  2. 


Boundaries 

Limiting  Statistics 

Corrections  to  the 
Limiting  Statistics 

Unchanged 

Unchanged 

Unchanged 

Replace  ♦/fi>l 

by 

Unchanged  except 
that  A is  not  known 
explicitly 

Replace  a by  a* 

Table  2.  Changes  Needed  to  Apply  the  Formulas  of 

Secs.(l)-(5)  to  the  Fully  Saturated  Regime. 
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The  higher  order  statistics  in  the  partially  saturated  regime 
are  more  complicated.  For  the  case  p < 2 everything  can  be  worked 
out  in  detail  and  the  results  will  be  given  in  the  next  section. 
However,  for  p > 2 the  path  integrals  yield  only  qualitative  infor- 
mation: It  is  summarized  in  Appendix  D. 

Finally,  in  multiple  scale  media  the  notion  of  multiple  Fermat 
paths  should  be  used  with  care.  They  exist  but  there  are  so  many 
of  them  that  they  cannot,  even  in  principle,  be  completely  resolved. 
Nevertheless,  the  notion  is  useful  in  interpreting  the  path  integral 
calculations  and  will  continue  to  be  employed. 


8 . The  Partially  Saturated  Regime  for  p < 2 


The  partially  saturated  regime  for  p < 2 is  of  considerable 
practical  importance.  Many  atmospheric  optics  experiments  lie  in 
this  region  and,  luckily,  the  complete  statistics  of  can  be 
worked  out.  There  is  a natural  small  parameter  g defined  by 


2-p 


(8.1) 


For  p = 5/3,  8 is  related  to  Tatarskii's  by 

8 = 1.19  and  to  the  intensity  fluctuations  as 

computed  in  the  Rytov  approximation  by/(£n  I)  £nlX  = 0.808 

Rytov 


-5/2 


The  signal  statistics  will  be  given  through  order  g. 

Partial  saturation  is  due  to  the  appearance  of  multiple  Fermat 
paths  all  lying  within  L of  each  other.  The  larger  scales  (~L) 
will  tend  to  correlate  the  locations  of  these  paths  leading  in 
general  to  a complicated  statistics.  However,  for  p < 2 the  spectrum 
is  so  heavily  weighted  toward  small  scales  that  the  locations  of 
the  Fermat  paths  turn  out  to  be  uncorrelated.  This  is  not  the 
case  for  p > 2 where  the  multiple  Fermat  paths  become  correlated 
and  the  path  integral  yields  only  qualitative  information  (see 
App.D) . Even  for  p < 2 where  the  locations  of  the  paths  are 
uncorrelated  the  large  scales  can  still  correlate  the  phases  along 
different  Fermat  paths.  We  will  see  this  at  the  end  of  the  section 
when  coherences  in  frequency  are  studied. 

Consider  Eqs.(3.1)  and  (3.2)  for  <I*>  in  the  partially 
saturated  regime  with  p < 2.  In  the  integration  region  (a)  the 
separation  between  members  of  a pair  of  paths  Vj^(2)  (using  the 
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notation  of  App.  B)  must  be  such  that  d(|vj^(z)|)  < 1,  i.e., 
|Vj^(z)|  < The  distance  ^2(2)  between  pairs  (again  in 

the  notation  of  App.B)  will  be  limited  by  the  oscillating  terms 
in  the  path  integral  to  values  such  that  n | V2  (z)  | | (z)  | ^ L or 

lv2(z)|  < Note  that  the  ratio  of  the  cutoff  on  IV2I  to 

that  on  |vj^|  is  and  is  large.  Now  both  and  IV2I  are 

small  compared  to  L and  Eq. (7.2)  can  be  used  to  evaluate  M in 
Eqs.(3.2)  or  (B.3).  Taking  account  of  the  fact  that  |Vj^|  <<  IV2I 
the  expression  for  M in  Eq.(B.3)  of  App. (B)  then  becomes^® 


2 -1 
M = $ R 


dz  - 4>^R 


^P(P-i)  f 1^ 


,2,-»-  ,p-2 

(z)  V.,  (z) 


dz  (8.2) 


and  when  |vj^/l|  ~ $ and  IV2/LI  ~ the  second  term  on  the 

right  hand  side  of  Eg. (8.2)  is  of  order  B and  can  be  dropped.  This 
is  the  same  thing  as  saying  that  different  pairs  of  paths  in  the 
path  integral,  or  equivalently  different  Fermat  paths  in  the  sense 
of  Sec. (6),  are  not  correlated  and  <I^>  becomes  2<I>^,  What  is 
happening  is  that  for  p < 2 the  fractional  power  behavior  of  d 
at  small  separations  is  making  the  arguments  of  Sec.  (3)  valid  even 
though  the  different  pairs  are  separated  by  less  than  L.  Note 
that  this  will  only  happen  for  p < 2.  The  generalization  to 
<l”>  is  straightforvard  and  the  result  is  a Rayleigh  distribution 
with  <l">  = n!  <I>". 

The  true  test  of  the  method  comes  when  one  evaluates  the 
corrections  to  Rayleigh  statistics.  It  is  shown  in  App.C  that  to 
order  B 


<l">  = nl<I>"(l  + I n(n-l)C(p)B) 


(8.3) 
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where  C(p)  is  a constant  which  depends  only  on  p.  This  constant 
is  evaluated  in  App. (B)  and  C(5/3)  = 1.06.  The  corrections  are 
small  for  small  B showing  that  the  approximation  scheme  is  con- 
sistent but  there  will  be  significant  deviations  from  a Rayleigh 
distribution  when  !/<!>  > /2/BC (p)  . 

The  statistics  of  as  a function  of  source  and 

receiver  locations  can  be  investigated  in  a similar  way.  One 
finds  that  they  are  Gaussian  and  at  equal  times  and  frequencies 
the  results  of  Secs.  (3)  and  (4)  hold  in  the  limit  3 = 0.  There 
are  coherence  tails  of  order  6.  These  are  discussed  in  App.C. 

In  the  fully  saturated  regime  the  dynamics  of  the  medium 
enters  only  through  D(t).  This  is  not  always  true  in  the  partially 
saturated  regime.  It  is  true  when  the  Taylor  hypothesis  is  valid 
(a  frozen  field  convected  by  a "wind")  and  the  statistics  in  time 
can  be  obtained  from  the  spatial  statistics.  However,  one  can 
consider  a different  kind  of  medium  where  the  time  dependence 
of  p is  associated  with  linear  waves  whose  dispersion  relation  is 
0)  ~ . The  Fourier  transform  of  the  second  time  derivative 

p of  p will  then  behave  like 

p(|q|)  = (const)  |q|  ” (8.4) 

at  large  |q|.  For  the  Taylor  hypothesis  Eq.(8.4)  holds  with 
6=2  and  in  general  6 can  be  considered  as  being  defined  by 
Eq.(8.4).  Assuming  p < 2,  the  statistics  of  £ at  unequal  times 
are  Gaussian  in  the  partially  saturated  regime  provided  that 
p - 6 <0.  This  can  be  verified  by  explicitly  computing  the 
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corrections.  For  p < 2 and  p - 6 < 0 the  corrections  to  Gaussian 
statistics  are  fractionally  small  for  small  S and  the  results  of 
Secs. (3)  and  (4)  continue  to  hold  at  unequal  times.  However  for 
p - 6 > 0,  a direct  calculation  shows  that  the  corrections  to 
Gaussian  statistics  are  not  fractionally  small  and  therefore  that 
the  approximation  scheme  of  Sec.  (3)  is  not  consistent  at  unequal 


times. 

To  see  what  is  happening  for  p - 6 > 0 one  can  compare  the  path 
integrals  for  <I(t')I(t)>  and  ^^£,*  (t ' ) ^^(t)  ^ . The  latter 

is  sensitive  to  the  time  dependence  of  the  phase  of  vl,  while  the  former 
is  not.  A rather  involved  but  straightforward  calculation  then 
shows  that  for  p - 6 > 0 the  signal  moves  more  rapidly  in  phase 
than  in  amplitude.  This  is  to  be  contrasted  with  the  case 
p - 6 < 0 where  the  time  statistics  are  Gaussian  and  according 
to  Eq.(1.18)  there  is  no  tendency  to  move  in  phase  as  opposed  to 
amplitude.  As  long  as  p < 2 the  signal  has  a Rayleigh  distribution 
and  over  a long  time  the  track  of  the  signal  will  fill  out  a disc 
in  the  complex  plane.  The  difference  between  p - 6 < 0 and 
p - 6 ->  0 comes  in  how  this  disc  is  filled  up.  For  p - 6 < 0 the 
signal  is  Gaussian  and  it  will  make  a track  of  the  type  shown  in 
Fig. (5a)  which  looks  something  like  a random  walk.  However, 
for  p - 6 > 0 the  track  will  wrap  around  in  phase  and  slowly  move 
in  and  out  in  amplitude  as  shown  in  Fig.  (5b). 

These  peculiar  features  of  time  statistics  in  the  partially 

saturated  regime  can  be  understood  in  terms  of  Fermat  paths.  We 

« i<t>.  (t) 

know  that  t(t)  is  schematically  EA,  (t)e  ^ where  the  locations 

k ^ 
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of  the  paths  are  uncorrelated  (for  p < 2)  but  the  large  scales 
may  correlate  the  phases  (j>j^(t).  The  question  of  random  walking 
vs  phase  wrapping  is  equivalent  to  the  question  of  whether  or 


not  the  time  derivatives  ^ (()j^ 


are  correlated.  For  p - 6 < 0, 


the  time  derivatives  are  sufficiently  weighted  towards  small  scales 
that  the  are  uncorrelated  and  the  signal  random  walks.  However 
for  p - 6 > 0,  the  effect  of  the  large  scales  is  strong 
enough  to  producv^:  a correlated  phase  derivative  common  to  all  the 
Fermat  paths. 

Propagation  of  sound  in  the  ocean  is  an  example  of  a situation 

7 

where  £(t)  phase  wraps  in  the  partially  saturated  regime.  For 

the  ocean  p = 2,  5 •>  0 and  in  this  special  case  it  is  possible  to 

work  out  the  detailed  statistics  of  (fct) . However,  for  other 

combinations  of  p and  6 it  is  not  possible  to  compute  fourth  and 

higher  moments  of  (£(t)  analytically,  except  when  p - 6 < 0. 

Checking  consistency,  it  was  stated  above  that  for 

p < 2 the  statistics  of  ^(r,r  ) as  a function  of  r and  r are 

o o 

Gaussian  in  the  partially  saturated  regime.  If  the  time  derivatives 

on  the  right  hand  side  of  Eq.(8.4)  were  replaced  by  spatial 

derivatives  we  would  have  6=2.  Since  p - 2 < 0 for  p < 2 it 

IS  consistent  that  the  statistics  in  r and  r are  Gaussian  and 

o 

that  the  statistics  in  time  are  Gaussian  when  the  Taylor 
hypothesis  (implying  6=2)  is  valid. 

For  p < 2 and  p - 6 < 0,  the  statistics  of  £ (r,r^,t)  in  the 
partially  saturated  regime  are  essentially  the  same  as  in  the 
fully  saturated  regime.  The  reader  may  therefore  wonder  what  the 
basic  distinction  between  the  regimes  is.  The  answer  turns  out  to 
lie  in  the  statistics  in  frequency. 
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Let  us  examine  the  path  integral  for  ^ ’ 

Up  to  a normalization  it  is 


<£*(u^)£(»2)£.*(a.3)e(»4)  > 


/ 


(paths)  exp 


R 

■4  .ui.  r 

2 j.i  = J J 


’dz  - N 


(8.5) 


where  with  the  Markov  approximation,  surpressing  time  t 

R 

. -v  -*■ 

dz 

■ i ] ■ 

0 


N = i L (-1)^"^^ 

i.j  = l 


K 

0) . Cl) . r ^ ^ 

- 2^-  J p ( I r^  (z) -r^  (z)  I ) ( 


(8.6) 


and  for  simplicity  the  medium  has  been  assumed  to  be  nondispersive . 
There  are  the  usual  two  important  regions  of  path  space  (a)  and  (b) . 
Let  us  concentrate  on  (a)  where  and  [r^-r^j  < . 

First  we  will  see  how  Gaussian  statistics  arise  in  the  fully 
saturated  regime  and  then  see  how  the  partially  saturated  case 
differs.  In  the  fully  saturated  case  typical  values  of,  say, 
jr^^-r^l  are  large  compared  to  L and  p(|r^-r^|)  can  be  set  equal 
to  zero.  Ignoring  correlations  between  the  different  pairs  then 
yields 

R 

Hz^l  I + E ) I p(|?.  (z)-?.  (z)  |)dz  (8.7) 

^ Vi,j  = l i,j  = 3/  J ^ ^ 

0 

which  is  a sum  of  two  terms  one  of  which  depends  on  and  0)2 

the  other  on  u),  and  co,  and  the  result  is  Gaussian  statistics.  In 
3 4 

the  partially  saturated  case  typical  values  of  Ir^^-r^l  are  small 
compared  to  L and  ^(Irj^-r^l)  is  approximately  equal  to  0(0).  Now 
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we  have  to  set  correlations  between  the  different  pairs  of  paths 
equal  to  p(0)  rather  than  zero  and  N becomes 


(8.8) 


The  path  integral  again  factors  into  a product  of  two  double  path 
integrals  and  is  expressable  in  terms  of  A as  defined  in  Eq.(2.22). 
Collecting  the  contribution  from  both  regions  (a)  and  (b)  and 
supplying  the  correct  normalization  yields 


exp  ^ E^(-l)^ 


((Jj^-U)^)  +A  A (0)^ 


(8.9) 


Because  of  the  common  exponential  factor  in  front  of  the  two  terms 

on  the  right  hand  side  this  is  not  Gaussian  statistics.  What  it 

corresponds  to  is  an  <£  of  the  form  ^(w) (u)  = e^“'^X(‘^)  where 

2 -2 

t);  is  a real  Gaussian  random  variable  with  <tp>  = 0,  <i()  > = R5(0)c 
and  x(<jj)  is  an  independent  complex  Gaussian  random  variable  with 
zero  mean  and  covariances  <^x  (t^)  X (<^ ' )^  = <^X*  (i^)  X*  (w ' ) ^ = 0 and 
^X*(w)X  (<*>')>  = A(  u)'-  (ij)  . It  is  straightforward  to  verify  that  this 


1 

I 
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ansatz  does  in  fact  yield  the  correct  2n-th  moment  of  (oj)  in 
the  partially  saturated  regime.  In  particular  the  second  moment 


<^exp 

(ui ' -m)  'y 

= 

<^exp 

[i  (u) ' -w)  (u))  X (ii>' ) y 

= 

exp 

(a)-a)')^R^(0)c"2j  A{m'-a))  (8.10) 

comes  out  right. 

For  a 

5 -2 

dispersive  medium  <ij;  > becomes  u)^  as  in 

Eq.  (1.  20)  and 

rather 

than  c appears  in  A. 

Thus  the  fundamental  distinction  between  the  fully  and 
partially  saturated  regimes  is  that  in  the  former  the  statistics 
in  frequency  are  Gaussian  while  in  the  latter  they  correspond  to 
a phase  times  a Gaussian.  Well  inside  the  partially  saturated 
regime  tOg  is  small  compared  to  the  width  in  oj  of  A,  The  phase 
eiwij;  dominates  the  moments  of  t except  for  correlations 

involving  only  ] (£  (w)  \ ^ where  4*  cancels.  As  the  boundary  = 1 

of  the  fully  saturated  regime  is  approached  the  width  of  A(oj) 
becomes  comparable  to  and  upon  passing  into  the  fully  saturated 
regime  A dominates  the  moments  and  the  signal  becomes  Gaussian. 

In  the  terminology  of  Sec.  (4),  for  partial  saturation  the  spread 
is  small  compared  to  the  wander.  In  pulse  propagation 
represents  a quasi-deterministic  wander  which  dominates  . 

The  phase  cancels  out  in  the  integral  (Eq.(4.8))  for  P(t) 

and  the  spreading  of  a pulse  is  proportional  to  the  inverse  width 
of  A. 
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, the  non- 


In  terms  of  Fermat  paths  C (lii)  = EA,  (a))e 

• k ^ 

Gaussing  statistics  can  be  understood  as  follows.  Each  <(>.  (w) 

IV 

can  be  written  as  cjtp  + where  uj\Ji  is  a common  phase  generated 

by  the  larger  scales.  The  phase  differencesA(J)  (to)  are  due  to 
the  small  scales.  They  vary  from  path  to  path  and  are  responsible 
for  the  Gaussian  factor  Note  that  only  correlations  in 

frequency  measure  directly.  Correlations  in  space  or  time  see 
only  ^ip  or  iji  which  for  p < 2 and  p - 6 < 0 are  dominated  by  small 
rather  than  large  scales,  leading  ultimately  to  Gaussian  statistics. 
The  phase  wrapping  in  time  for  p - 6 > 0 is  a remnent  of  \p. 

The  statistics  of  ^ in  the  partially  saturated  regime  are 
summarized  in  Table  3. 


1 

Intensity 

Distribution  j 

Variations  in 

space 

1 time 

L.  J 

frequency 

Rayleigh 

1 

1 

Gaussian 

p - 6 < 0 

p - 6 > 0 

phase  times 
a Gaussian 

Gaussian 

1 

phase 

wrapping 

1 

Table  3.  The  Statistics  of  In  the  Partially 
Saturated  Regime  for  p < 2. 


The  reader  may  be  curious  as  to  what  happens  at  p = 2.  The 

"small"  parameter  6 is  then  equal  to  unity  but  according  to  App.  C 

the  coefficient  C(2)  in  Eq.  (8.3)  vanishes.  A detailed 
7 

investigation  then  shows  that  the  corrections  to  Rayleigh 

statistics  in  the  partially  saturated  regime  are  of  order  (iln*)”^. 

More  generally,  if  p = 2 and  ini>  is  large  the  statistics  given  in 

Table  3 apply  with  errors  of  order  (iln<!>)~^.  At  p = 2 it  is 

possible  to  compute  A.  It  is  given  by  Eq.(1.21)  with  replaced 
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r 

-i 

I 

1 

! 

by  In  general,  a medium  with  |p-2|£n<I>  < 1 will  act 

like  one  with  p = 2. 

f As  mentioned  before  the  case  of  partial  saturation  for 

p > 2 is  discussed  in  App.  D. 


9.  INHOMOGENEOUS  AND  ANISOTROPIC  MEDIA 
In  practice,  random  media  are  only  locally  homogeneous  and  the 
covariance 

p(x  - x',  t - t‘;  x)  = < p(x  ,t)  p(  x',  t' ) > - < fx(  X , t)  > < (i(x  t')  > (9. 1) 

depends  on  position  x = ^ {x  + x' ) . It  is  always  assumed  that  the  variations 
of  p in  x-x'  are  much  more  rapid  than  those  in  x but  over  a long 
propagation  path  the  dependence  on  x cannot  always  be  neglected.  Also, 
in  an  inhomogeneous  medium  < |jl(x)  > = Pq(x)  will  generally  not  be  a 
constant  and  consequently  cannot  be  absorbed  in  the  definition  c = £o/k  . 
Finally,  the  medium  can  be  statistically  anisotropic  so  that  p depends 
on  the  orientation  of  x-x'  as  well  as  its  magnitude. 

To  obtain  tractable  path  integrals  in  an  inhomogeneous  medium 
we  will  have  to  approximate  the  path  dependence  of  x in  p by  evaluating 
X along  some  central  path  which  will  turn  out  to  be  an  unperturbed  ray. 
From  Sec.  (3)  we  know  that  paths  are  separated  by  L (the  precise 
definitions  of  L,  <I>  and  for  inhomogeneous  anisotropic  media  will  be 
given  below)  and  the  problem  will  be  tractable  if 

(i)  for  changes  in  x of  order  L$/D  the  corresponding  variations 
in  p can  be  neglected.  It  will  also  turn  out  to  be  necessary  to  expand 

in  powers  of  distances  between  paths  and  we  will  have  to  require  that 

(ii)  M'q(x)  is  slowly  varying  over  distances  of  order  LiJ/D  . 

The  one  other  condition  is  that 

(iii) the  parabolic  wave  equation  is  a valid  approximation. 
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When  M’q(x)  is  not  a constant  this  requires  .that  the  normals  to  the  wave 
fronts  in  the  unperturbed  problem”  where  p.(x)  = ^.^(x)  remain  close 
to  the  z-axis.  If  this  is  true  locally  but  not  globally,  then  solutions  based 
on  the  parabolic  approximation  can  be  patched  together  in  the  obvious  way. 
When  conditions  (i) , (ii)  and  (iii)  are  met  it  is  reasonably  straight- 
forward to  extend  the  path  integral  method  to  inhomogeneous  and  aniso- 
tropic media.  It  amounts  to;  (1)  showing  that  with  suitable  definitions 
of  3.  and  D,  <C  > remains  exp[-i$^]  and  Eq.  (1.14)  continues  to 
hold,  (2)  finding  a suitable  definition  for  f2  and  then  showing  that  the 
boundaries  of  the  saturated  regimes  are  still  given  by  Table  1,  (3)  showing 
that  in  the  fully  saturated  regime  the  statistics  of  ^ are  Gaussian  and 
that  in  the  partially  saturated  regime  they  are  (for  p < 2)  as  given  in 
Table  3,  (4)  giving  new  formulas  for  the  corrections  to  Gaussian  statistics 
and  coherence  tails  and  (5)  giving  a method  for  computing  A(u)).  These 
steps  will  be  carried  out  in  order.  In  doing  so  it  will  be  assumed  that  a 
ray  approximation  is  valid  for  the  unperturbed  problem  with  |i.  = 

A.  The  First  and  Second  Moments 

The  path  integral  for  < C > will  contain  a factor 


r ,^2  R R 

exp  ~ Y f I Pr  -r(z')  + e (z 
'-  0 0 ' ^ 


z')  .0;  i(r(z)  + r(z'))  + |e  (z+  z')j  j 

(9.  2) 


The  path  dependence  of  the  third  argument  of  p will  be  approximated  by 
setting  j{r{z)  + r(z')  ) = 3(z)  where  z = 5(z+  z')  and  s is  the  unperturbed 
ray  satisfying 
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8"(z)  + V (A  (s(z)  + e z)  = 0 
u z 


(9.3) 


Here  s = (s  , s ) is  a two  dimensional  vector  and  (s(0),  0)  and  (s(R),  R) 
X y 

are  the  source  and  receiver  co-ordinates.  If  there  is  more  than  one  un- 
perturbed ray  connecting  the  source  to  the  receiver  it  is  assumed  that 

they  are  far  enough  apart  that  the  path  integral  reduces  to  a sum  of 

PI 

(statistically)  independent  terms  coming  from  paths  near  each  ray. 

Defining  a new  path  u(z)  by  r(z)  = s(z)+u(z),  the  Markov  approxima- 
tion now  amounts  to  setting 


r(z)-r(z')  + e (z-z')  ~ (s'(z)+e  )(z-z') 

z z 


(9.4) 


The  essence  of  the  approximation  is  neglecting  u(z)  - u{z').  By  require- 

ment  (ii)  the  substitution  s{z)  - s(z' ) st  s ' (z)  (z  - z' ) is  always  valid. 

The  reader  will  note  that  by  (iii)  s'  is  actually  small  compared  to  e . 

z 

However,  in  a sufficiently  anisotropic  medium  s'  cannot  be  neglected 
on  the  right-hand  side  of  Eq.  (9.  4).  Assuming  for  the  moment  the  validity 
of  the  Markov  approximation,  the  analog  of  p(0,  0)  will  be  p(0, 0;  z) 
where 

oo 

p(0 , 0 ; z)  = f p(  (s'(z)  + e ) u,  0;  s(z)  + e z)  du  (9.  5) 

.CO  ^ ^ 

and  the  path  integral  for  <^>,  which  is  now  trivial  since  p no  longer 
contains  the  path  u,  will  yield  <^>  - exp  [ - j where 

- R 


« = k'^  J p(0,0;z)dz 

0 


(9.6) 
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Continuing  to  assume  the  validity  of  the  Markov  approximation 
the  next  thing  to  compute  is  *(2)  S (1)  >,  There  are  two  paths 
r^  = s + Wj  and  ^2  " ® ^2  where  s satisfies  Eq.  (9.  3)  with  the 

boundary  conditions  s(0)  = 2(^01  "^^02^  ^ and  the 

approximation  is 
R R 

Jdz  J dz'p(7.(z)  -T  (z')  + e (z-z').  t -t.;i(T.(z)  +T.(z'))+Ae  (z+z')) 

0 0 ■’  1 J 1 J z 

R 

« f p(w.(z)  - W.(z);t.  - t.:z)dz  (9 

^ 1 J 1 J 

for  i,  j = 1,  2 where 


(9.7) 


p(w,  t;  z)  - r p(w+(s'(z)  +e  )u,t;s(z)  +e  z)du 

7.  ’7 


(9.8) 


The  path  integral  for  <<^*(2)  (f(l)  > is  then 


<^*(2)^(1)>  = J d^(paths)  expTi  S (path  1)  -i  S (path  2) 
4k  L 


ix  1 

- rd(w  (z)  - w (z),  t -t,;  z)dz 
J I 2 12  J 


where 


0 " (^'(2)  )^  - |iQ(r(z)  + e^z  )ldz 


(9.  9) 


(9.10) 


d(w,  t;z)  = k [ p (0.  0;  z)  - p(w.  t;  z)  ] 


(9.  11) 


Introducing  paths  u = 2(Wj+w2)  and  v = we  can,  according  to 
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(ii),  expand  SQ(pathl)  -SQ(path2)  in  powers  of  u and  v and  keep 
only  the  leading  terms  which  are  quadratic.  Proceeding  in  this  way 


yields 

I* 


22 


<r^(2)<f(l) — _ |2k^  I [d^(paths)expr ik[(u'(z)  • v'(z)  - u.(z)  v.(z)  p..(z)  ) dz 

^ I-  i ^ J 


R 

-f  d(v(z),  z)dz  J 


(9.12) 


where  the  two-by-two  matrix  (in  e - e space)  u (z)  is 

X y ij 


ij  dx.  9x,  0 * ^ » 

1 j x=s(z)+ez 

z 


(9.13) 


The  path  u now  appears  only  as  a linear  factor  in  the  exponential  and 
integrating  over  it  will  produce  a product  of  6-functions  which  force 
v(z)  to  be  equal  to  the  special  path  z^z)  which  satisfies  the  differential 


equation  and  boundary  conditions 
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V'."(z)  +p,.(z)t>',  (z)  = 0 

1 J 


T?(0)  = r -7 
' ' 01  02 


t.(R)  = 


(9. 14) 


Then  setting  v equal  to  ij'  in  d the  remaining  path  integral  just 
17  Z 

produces  (2k  ^ q\  and  one  finds  Eq.  (1.  14) 

<<f*(2)^(l)>  = Qz)  tf^d)  exp[-iD]  (1.14") 


with 


R 


D=  Zj  d(x^(z) , t,  - t,  ; z)dz 
0 * ^ 


(9. 15) 
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The  object  D defined  in  Eq.  (9.15)  is  just  the  phase  structure 
function  of  first  order  geometric  optics^' ^ for  a general  inhomogeneous 
anisotropic  medium  which  satisfies  (i) , (ii)  and  (iii) . Note  that 
is  always  linear  in  - r^^  and  When  is  a constant, 

= (r^^  - (R  - z)  /R  + (r^  "^2^  ^ homogeneous  iso- 

tropic medium  Eq.  (9.  15)  reduces  to  Eq.  (1.  16). 

For  an  isotropic  medium  where  p depends  only  on  the  magnitude 
of  x-x'  the  Markov  approximation  is  valid  whenever  the  parabolic 
wave  equation  is.  The  reason  is  the  same  as  in  Sec.  (Z).  In  App,  (E) 
the  formula  for  the  first  correction  to  the  Markov  approximation  to 
<&*t  > is  given.  One  can  explicitly  verify  that  the  error  is  small 
when  the  parabolic  wave  equation  is  valid. 

The  situation  for  anisotropic  media  is  more  complicated.  Consider 
an  anisotropic  but  homogeneous  medium  with  constant  p^.  Typical  in- 
homogeneities will  not  be  spherically  symmetric  and  one  needs  to 
consider  the  three  cases  shown  in  Figs.  (6a),  (6b)  and  (6c).  The 
asymmetric  i n h om  o g e n e i t i e s introduce  a new  small 
angle  0^,  the  ratio  of  the  small  dimension  to 
the  large  one.  Examining  the  error  In  the  Markov  approximation  as 
given  in  App.  (E)  one  finds  that, for  the  case  shown  in  Fig.  (6a),  the 
Markov  approximation  fails  when  the  r.  m.  s.  multiple  scattering  angle 
is  of  order  . For  the  case  shown  in  Fig.  (6b),  it  fails  when  the 
r.  m.  s.  multiple  scattering  angle  is  of  order  of  the  angle  of  incidence 
0^  and  for  the  situation  in  Fig.  (6c),  it  fails  when  the  r.  m.  s.  multiple 
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scattering  angle  is  of  order  unity,  i.  e.  , 


when  the  parabolic  wave  equa- 


tion fails.  Since  0^  can  be  small  compared  to  unity  the  Markov  ap- 
proximation can  fail  in  an  anisotropic  medium  before  the  parabolic 
wave  equation  does  but  only  for  some  propagation  paths.  When  it  fails 
Eq.  (1. 14)  is  not  valid  and  this  represents  a defect  in  the  theory  which 
is  not  easy  to  remove. 

It  should  not  be  surprising  that  the  Markov  approximation  can 
fail  sooner  in  an  anisotropic  medium.  The  Markov  approximation  can 
be  interpreted  as  the  statement  that  the  system  has  "no  memory"  in 
range,  i.  e.  , that  scatterings  at  a given  range  point  are  independent  of 
previous  distant  scatterings.  In  an  isotropic  medium  this  will  be  true 
as  long  as  the  r.  m.  s.  multiple  scattering  angle  is  small  and  the  wave 
keeps  moving  in  the  same  direction.  However,  in  an  anisotropic  medium, 
when  the  scattering  by  a given  inhomogeneity  can  be  highly  dependent  on 
the  angle  of  incidence,  a distant  scattering  which  has  deflected  the  wave  only 
through  a small  angle  will  not  be  "forgotten."  For  the  inhomogeneities 
shown  in  Fig.  ( 6 ) the  scattering  is  strongly  dependent  on  angle  of 
incidence  (measured  from  the  long  axis  of  the  inhomogeneities)  when 
the  angle  is  of  order  0^.  When  the  incident  wave  is  along  the  long  axis 
as  in  Fig.  (6a),  it  begins  to  remember  previous  scatterings  when  the 
scattering  angle  builds  up  to  0^  and  the  pieces  of  the  wave  have  inci- 
dence angles  greater  than  0^.  For  the  case  shown  in  Fig.  (6b)  the  past 
history  of  the  wave  becomes  important  when  pieces  of  the  wave  have 
been  deflected  by  0^  and  are  incident  along  the  long  axis.  When  0^ 
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approaches  ir/2  as  in  Fig.  (6c)  the  process  has  no  memory  as  long 
as  the  r.  m.  s.  multiple  scattering  angle  i«  less  than  unity. 

Yet  another  way  to  understand  the  peculiarities  of  anisotropic 
media  is  to  return  to  the  remarks  following  Eqs.  (6.  8)  and  (6.  9) . For 
an  isotropic  medium  the  average  of  |1  integrated  along  a Fermat  path 
w will  be  the  same  as  the  average  of  p integrated  along  the  unperturbed 
ray  s as  long  as  the  r.  rn.  s.  multiple  scattering  angle  is  small.  How- 
ever, in  an  anisotropic  medium  the  average  of  p integrated  along  a 
path  can  be  very  sensitive  to  the  local  direction  w'  of  the  path.  In 
fact,  for  the  situation  shown  in  Fig.  (6a),  the  average  of  p integrated 
along  a Fermat  path  deviates  from  the  average  along  an  unperturbed  ray 
as  soon  as  1 ^ | ~ 6^  and  for  the  situation  in  Fig.  (6b)  when  ] w ' | ~ 9^, 
This  leads  to  the  same  criteria  as  before. 

The  combination  of  an  anisotropic  medium  and  a spacially  varying 

leads  to  a new  set  of  complications.  This  will  be  illustrated  for 
propagation  in  a channel  where  the  unperturbed  rays  make  loops  as  shown 
in  Fig.  (7)  and  where  the  long  axis  of  the  inhomogeneities  is  parallel  to 
the  channel  axis.  The  medium  will  also  be  assumed  to  be  statistically 
homogeneous  in  the  direction  of  the  channel  axis  but  not  necessarily  in 
the  transverse  directions.  (This  is  a prototype  of  the  physical  situation 
which  Occurs  for  sound  propagation  in  the  ocean.  ) The  scattering  will 

be  strongest  when  the  tangent  to  the  unperturbed  ray  is  pointing  along 
the  long  axis  of  the  inhomogeneities,  i.  e.  , at  the  turning  points. 
For  small  0^  one  can  in  fact  ignore  all  of  the  propagation  path  except 
for  a set  of  discrete  regions  around  turning  points  where  the  tangent  to 
the  ray  is  within  0^  of  the  channel  axis.  Assuming  that  a Markov 
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approximation  is  valid  for  propagation  through  one  of  these  regions,  it 

will  also  be  valid  for  propagation  through  many  turning  points  provided 

only  that  the  average  scattering  at  a given  turning  point  is  at  most 

weakly  dependent  on  scatterings  at  previous  turning  points.  Assuming 

that  the  turning  points  are  separated  by  more  than  a coherence  length 

the  effect  of  previous  scatterings  will  be  a random  modulation  of  the 

range  and  (transverse)  location  s ^ of  a turning  point.  Now  the 

average  scattering  around  a turning  point  is  dependent  only  on  its 

location  s^  in  the  channel  and  not  on  its  range  z^.  Thus  the  Markov 

approximation  will  be  valid  out  to  range  such  that  random  variations 

in  Sq  are  big  enough  to  change  the  average  scattering.  This  turns 

7 

out  to  be  a much  longer  range  than  that  for  which  the  r.  m.  s.  multiple 
scattering  angle  (which  is  dominated  by  variations  in  z^)  becomes  of 
order  6^.  The  extended  validity  of  the  Markov  approximation  can  be 
demonstrated  explicitly  using  the  Fermat  path  formalism  of  Sec.  (6). 
One  works  out  the  properties  of  Fermat  paths  which  are  randomly  de- 
flected at  turning  points  and  then  compares  averages  of  p integrated 
along  these  paths  to  averages  of  p integrated  along  the  unperturbed 
ray.  For  a given  channel  one  can  then  find  out  when  the  Markov  ap- 
proximation will  break  down.  The  result  is  just  the  criteria  stated 
above. 

B.  The  Saturated  Regimes 

It  will  temporarily  be  assumed  that  the  medium  has  a single 
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scale.  Then  in  an  anisotropic  Inhomogeneous  medium  the  scale  length 

t L becomes  a z-dependent  two-by-two  matrix  (in  e^  - e space)  defined 

..  22 

by  the  expansion  of  p 

p(w,0:z)  = p(0,  0;  z)  [ 1 - I (L“^{z)  ).^  w^  w^  + 0([  wp)  ] (9.16) 

The  first  task  in  discussing  the  saturated  regimes  is  to  find  the 
correct  definition  of  and  establish  their  boundaries.  The  general 
definition  of  n will  involve  L and  some  geometric  parameters  asso- 
ciated with  the  unperturbed  problem.  From  Secs.  (3)  and  (7)  one  can 
see  that  Si  measures  the  rate  at  which  the  phase  of  the  oscillating 
factor  in  the  path  integral  varies  as  a path  moves  away  from  an  unper- 
turbed ray.  To  examine  this  in  more  detail  consider  paths  that  leave 
the  source  at  z = 0,  go  to  the  receiver  at  z = R and  at  some  point  z^ 
in  between  are  separated  from  the  unperturbed  ray  s(  z^)  by  I . Let 
Sl^.z^)  be  the  minimum  of  SQ(path)  - SQ(unperturbed  ray)  taken  over 
all  paths  of  this  class.  The  minimum  is  achieved  for  a path  that  follows 
an  unperturbed  ray  from  the  source  to  (s(Zq)  + t,  z^)  and  then  another 
unperturbed  ray  from  (s(z^)  + f , z^)  to  the  receiver.  When  is 
a constant  S (f  , z^)  is  simply 

♦ = id*)^  B(Zq)  (9.17) 

where  B which  already  appeared  in  Eq.  (6.  11)  is 


-1  -1  -2 
and  S2  is  the  average  of  B L , i.  e. 


S2~^  = L~^  B 


(9. 19) 


Thus  fl  is  a measure  of  the  phase  change  required  to  move  a path  a 
distance  L away  from  the  unperturbed  ray.  In  general  there  is  a 


two-by-two  matrix  B defined  by  the  expansion  for  small  f 


•-22 


(9.  20) 


and  ^ will  be  an  average  of  L ^ B \ It  is  convenient  to  weight  the 


average  by  p(0,  0;  z)  and  fl  will  be  defined  as 

^ .2 


22 


= 7 (L^)..  (B'S..  = ^ 

2 13  V3  2 


f p(0, 0;  z)  (L  (z)  )..  (B'\z)  )..  dz 

1 0 


V - 

J p (0  , 0;  z)  dz 


(9.  21) 


With  this  definition  of  f2  one  can  follow  through  the  arguments  of  Sec.  (3) 

and  verify  that  saturation  and  Gaussian  statistics  are  expected  when 

$ > 1 and  ^/i2  >1.  A more  precise  procedure  is  to  compute 
3 3 3 

( <1  > - <1>  )/<!>  in  the  Rytov  approximation  to  find  the  boundary 

of  the  saturated  regime  and  then  in  the  saturated  regime  compute  the 

corrections  to  Gaussian  statistics  and  verify  that  they  are  small.  A 

straightforward  evaluation  of  ( < I^>  - <I>^)  /<  I shows  that  it 

' Rytov 

does  in  fact  exceed  unity  when  $ > 1 and  5>/Sl  > 1 indicating  that  the 

boundary  is  correct.  Using  the  formulas  of  App.  F one  can  verify  that  in 
the  saturated  regime  the  corrections  to  Gaussian  statistics  are  indeed 
small. 
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With  the  appropriate  change  in  the  definition  of  {L  )_  for 
p < 2,  the  same  procedure  can  be  extended  to  media  with  multiple 
scales.  The  result  is  that  with  f2  defined  as  in  Eq.  (9.  21)  the 
boundaries  given  in  Table  1 remain  correct  and  that  for  p 2 the 
partially  saturated  statistics  given  in  Table  3 also  remain  correct. 
To  actually  calculate  B_(z)  the  following  result  is  useful. 

pp 

Define  a Green's  function  g^A^.z  ) by 

o2 

— ^ g (z,  z’)  + p.j^(z)  z')  = 6_6(z-z') 

dz 

g^j(0,  z')  = g.j(R.  z')  = 0 


Then  it  is  straightforward  to  verify  that 


(B  (z)).,  = -g..(z,z) 

ij  ij 


C.  Correlations  in  Frequenci 


In  general  one  can  write 


— — = exp  --  A(cj-w') 


(1.19') 


(9.  22) 


(9.  23) 


where  the  exponential  factor  comes  from  geometric  optics  and  A is 
to  be  computed  from  the  path  integral.  The  geometric  optics  decor- 
relation frequency  w is 

u ^ \ f ^ r (z)  + e z,  t ; - (s  (z)  + e z,  t ) >1  dz^  \ (9.  2 

g L ^ ^ * w w z j y / 

and  for  a general  dispersive  medium  both  p and  the  unperturbed  ray  s 
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will  depend  on  w.  For  a nondispersive  medium  w is  equal  to  w/'I> 

where  $ is  evaluated  at  the  central  frequency  w. 

As  before  the  path  integral  for  A is  tractable  only  for 

media  with  p > 2.  The  derivation  or oceeds  as  in  Sec.  (3)  and 

^ / k c \\ 


after  introducing  scaled  paths  ^ = ^2(cj  - 3)  j ^ where 
(U  = and  c is  the  group  velocity  at  w = w the  path  integral 

for  A gives 


A(w)  = 


K(cj) 

K(0) 


(9.  25) 


where 


r R ^ 2 T 

K(a))  = Jd(paths)  exp  -i  J [(|''(z))  - ^^(z)  (|j._(z)  + i w h_(z)  )]  dzj 


(9.  26) 


with 


h..(z)  = c"^  p(0,  0;  z)(l"^(z)  ) (9.27) 

ij  g 

If  the  path  integral  for  K is  written  out  in  its  finite  form  it 
becomes  an  ordinary  integral  of  large  dimension  whose  integrand  is 
the  exponential  of  a quadratic  form.  Such  an  integral  is  proportional 
to  one  over  the  square  root  of  the  determinant  of  the  quadratic  form 
and  in  particular  A will  be  the  square  root  of  the  ratio  of  two  determi- 
nants. As  the  number  of  integration  points  goes  to  infinity  the  deter- 
minants become  functional  determinants.  There  are  two  equivalent 
methods^  for  computing  the  ratio  of  these  functional  determinants. 

In  the  first  method  one  has  to  find  all  the  eigenvalues  of  the 


differential  equation 
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0 


(9.  28) 


'^{n)  '*■(11) 

subject  to  the  boundary  conditions  ^ (0)  = | (R)  = 0.  Having  done 

this  A(w)  is 


A{oj)  = [ "JT  i- 

\ n , . ; 


1 + 1 — 
bJ 

n 


(9.  29) 


In  the  second  method  one  defines  a two-by-two  matrix  M..{z,  w)  by 
the  differential  equation 


■Ml'.(z,w)  + ^jL  (z)  M,  .(z,w)  + ibjh.,  (z)  M,  .(*■<*>)  - 0 
ij  ik  kj  ik  kj 


(9.30) 


and  boundary  conditions 


in  z 


M,.(0,  w)  = 0 


M!.(0,«)  = 6 

iJ  ij 


(9.  31) 


Then  A is  given  by  the  ratio  of  determinants 

^det  M(R,  0) 


AM  = 


det  M(R,  w) 


(9.32) 


As  an  example  of  how  A is  computed  consider  a homogeneous 

isotropic  medium  where  p..  = 0 and  h.  = 6 c'^p(0,  0)l'^  The 

H ij  g 

eigenfunctions  of  the  operator  in  Eq.  (9.  28)  are  then  of  the  form 


6 sin(n  Tr  z/R)  and  6.^  sin(n  tt  z/R)  and  the  eigenvalues  are  -n  tr  cj 
2 2 2 2 - 

and  -n^  TT  where  = c^L  /R  p(0,  0)  as  in  Sec.  (4).  The  infinite 
product  in  Eq,  (9.  29)  is  then  a product  over  two  sets  of  integers 
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A(oj) 


00 


00 


( TT  ^ 

— TT 

Vn-l  1 i - 

1 1 

V r ^ ^ ‘ ? 2 

■ ^^2= 

Hi  IT 


(9.  33) 


and  the  two  equal  factors  just  cancel  the  square  root.  The  result  is 


A{^)  = I 


= TT 


s.nM 

n n n ^ 0 \ nir  \ w^/ 


(9.  34) 


and  with  u)^  = Eq.  (1.  21)  is  reproduced.  To  compute  A the 


second  method  one  finds  immediately  that  M..(z,cj) 

ij 


(W  \2  / / . \-l\ 

1 \ . / z f iw  y \ 

— smi  - — ] , 

lu)/  yR\Wj^/  j 


M..(z,  0)  = 6.,z  and  Eq.  (9.32)  yields  the  expected  answer, 
ij  ij 

Once  CO  and  A have  been  determined  everything  proceeds  as 
g 

o 


in  the  homogeneous  isotropic  case.  In  particular  C satisfies  Gaussian 
statistics  in  the  fully  saturated  regime  and  in  the  partially  saturated 
regime  for  p < 2 (where  A is  unfortunately  not  known)  it  is  a phase 
times  a Gaussian.  There  is  one  new  point  worth  mentioning.  In  the 
calculation  of  Sec.  (4)  < "^(t)  > vanished  for  t < 0 because  A was 

analytic  in  the  upper  half  plane.  When  M--j(z)  is  non-zero  there  can 

be  a finite  number  of  positive  eigenvalues  co^ . Then  A is  no  longer 

analytic  in  the  upper  half  plane  and  <'^  (t)  > is  nonvanishing  for  t < 0. 

7 

This  in  fact  happens  for  propagation  of  sound  in  the  ocean. 

The  transition  to  inhomogeneous  anisotropic  media  has  now  been 

completed.  The  reader  who  is  interested  in  seeing  how  the  method 

works  in  detail  for  a realistic  problem  can  consult  the  book  of  Flatte, 

, 7 

et  al. 


64 


10. 


Conclusions 


The  path  integral  has  turned  out  to  be  a powerful  tool. 
It  has  provided  a precise,  (very  nearly)  complete  and  global 
picture  of  what  goes  on  in  the  saturated  regimes.  The  unsat- 
urated regime  where  the  Rytov  approximation  is  valid  could 
also  be  treated  by  path  integral  methods.  While  this  would 
lead  to  a more  unified  picture,  in  the  end  it  would  only 
amount  to  a rederivation  of  the  Rytov  approximation.  A more 
1'. uitful  endeavor  would  be  to  make  an  attack  on  the  remaining 
unsolved  problems  in  the  saturated  regimes.  For  situations 
where  a scalar  wave  equation  is  sufficient  and  the  (multiple) 
scattering  angles  are  small  the  remaining  problems  are: 

a)  How  to  compute  (except  numerically)  the  coherence  in  fre- 
quency, A(o),  for  multiple  scale  media  with  p < 2. 

b)  What  are  the  detailed  (beyond  those  given  in  App.(D)) 
statistics  of  ^ in  the  partially  saturated  regime  for 
4 > p > 2? 

c)  How  to  compute  the  second  moment  (^^*{2)  €{1)^  for  those 
propagation  paths  in  highly  anisotropic  media  where  the 
Markov  approximation  is  not  valid? 

d)  Wiiat  is  the  detailed  behavior  of  ^ at  the  boundaries 
between  the  unsaturated  and  saturated  regimes  and  between 
the  fully  and  partially  saturated  regimes? 

These  are  difficult  problems  which  may  not  have  any 
simple  solution  and,  in  particular,  the  path  integral  may  not 
be  the  best  method  for  attacking  them.  On  the  other  hand,  it 
is  quite  remarkable  that  the  use  of  Feynmann's  path  integral 
has  reduced  the  problem  to  a few  unknowns  which  occur  only  in 
special  cases. 


Among  the  other  methods  for  treating  wave  propagation 
in  random  media,  the  most  powerful  ones  use  the  Markov  ap- 
proximation from  the  beginning.  With  the  Markov  approximation 
one  can  derive  local  partial  differential  equations  for  the 

moments  of  c . ” These  equations  have  been  studied  exten- 

3 4 

sively,  especially  by  the  Russian  school.  ’ In  the  Markov  ap- 
proximation the  path  integrals  for  the  moments  are  formal  solu- 
tions to  these  partial  differential  equations.  The  equations 
for  the  first  and  second  moments  can  be  integrated  analytically 
and  correspondingly  the  path  integrals  can  be  done  analytically. 
For  the  higher  moments,  the  differential  equations  have  yielded 
only  some  information  about^  ^I(l)  1(2)^  . The  reason  that 
this  approach  has  not  yielded  more  is  that  to  determine  the 
asymptotic  (long-range)  behavior  of  a function  from  its  defin- 
ing partial  differential  equation  is  highly  non-trivial.  The 
path  integral  has  the  advantage  that  it  works  on  a global  rather 

than  local  level,  making  it  easier  to  determine  the  asymptotics. 

12 

The  reader  who  is  familiar  with  Mercier's  treatment  of 
the  phase  screen  problem  (an  idealized  case  where  all  the  scat- 
tering takes  place  on  a thin  sheet)  will  have  noticed  the 
similarity  between  his  methods  and  those  used  here.  The  simi- 
larity is  partly  just  the  mathematics  of  manipulating  integrals 
but  there  is  also  a physical  reason.  Any  medium  can  be  approxi- 
mated by  a (perhaps)  large  but  finite  number  of  phase  screens. 
The  wave  field  can  then  be  expressed  as  a large  but  finite 
dimensional  integral  over  the  surfaces  of  the  screens.  But 
this  is  just  the  path  integral  in  its  finite  form.  Thus  the 
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I path  integral  can  be  thought  of  as  a scheme  where  one  approxi- 

mates the  medium  by  n phase  screens  and  then  letting  n go 
to  infinity  recovers  the  original  problem. 

f 
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APPENDIX  A CORRECTIONS  TO  "niE  MARKOV  APPROXIMATION 


The  exact  path  integral  for  ((S*(2)5(l))  can  be  expanded  as 

<(J*(2)5(1))  = ^ ^ ^ /d^( paths)  exp 

m=0  0 L ' 


X [V  - V] 
o 


dz  - V 


(A.l) 


where  V is  given  by  Eq.  (2.10)  and  V is  the  Markov  approximation 

o 

given  by  Eq.  (2.11).  The  m = 0 term  is  the  Markov  approximation  and 
the  m = 1 terra  will  be  computed  below  for  the  special  case  ((S*(l)(!(l)) . 

For  (5*(1)(S(1))  = (l)  the  first  correction  contains  [v^  - v] 

which  can  be  replaced  by  -V  because,  as  may  be  seen  from  Sec.  (2) 

the  piece  proportional  to  V vanishes.  We  then  have  to  first  order 

o 

R 

(l)  = d^(paths)  exp  ik^v'(z)  • u'(z)dz  “ V (A. 2) 

0 1+k  I 0 

where  paths  v(z)  = r^(z)  - ^^(z)  and  u(z)  = J|r^(z)  + have  been 

introduced.  In  terms  of  the  Fourier  transform  p of  h (Eq.  (5.3))V  is 
R R 

aC  / A ~uR 2'  \ ■ "2' 

V =2k^dz  /dz  /d  qdq  P\\q  + q ,0)  e 
Q 1 0 ^ z z 

iq  • (u(z^)  - 

X e sin  Uq  • v(z^)  sin U q • v(z2)  (A. 3) 

where  q = (q  .q  ) is  a two-dimensional  vector.  Writing 
X y 


iq  • ( u(z  ) - u(z  ))  f -»  /-» 

» '1  2 ' = exp  iq  • / u 


(z)  [6(z  - - 6(z  - z^)']dz 
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and  inserting  V as  given  by  Equation  (A. 3)  into  Eq.  (A. 2) 

one  finds  that  since  V depends  only  on  v the  integral  over  u(z)  can 

o ^ 

be  done  and  that  it  leads  to  a 6-  functional  which  forces  v to  satisfy 
the  equation 


v^'(z)  = ^[6(z  - z ) - 6(z  - z„)] 


(A.4) 


with  the  boundary  conditions  v(0)  = v(R)  = 0.  In  terms  of  the  Greens 

function  g(z,z')  which  satisfies  — g(z,z')  = 6 (z  - z')  and  g(0,z^)  = 

/ ->  d z 

g(R, z ) = 0,  v(z)  is  constrained  to  be  v (z)  where 

o 


- e(z.Z2^) 


(A.  5) 


The  path  integral  is  then  done  by  replacing  v(z)  by  v (z)  in  V and  in 

o o 

the  representation  of  V and  the  final  result  is 


<«  - 


R R 


^/dz  /dz  /d 
0 Iq  Z*' 


2 , '%^"i  - ^2^ 

qdq  P (V  q + q ,0)  e 
z z 


X sin 


|jq  • v^(z^)j  sinjiq  • exp  v^(z)  | , oj  ^ 

^ r\ 


(A.  6) 


To  estimate  the  size  of  the  integral  in  Eq.  (A. 6),  one  notes  that 

|z  - z I will  be  restricted  to  order  L or  less  and  that  for  |z  - z ] ~ L, 
^ 2 12 
• v^z)  is  of  order  q L/k  q/k  which  is  assumed  to  be  small.  The 

sines  can  then  be  expanded  and  using  the  identity 


[g(Zl,Zl)  - s(Z‘^,z^^'>'][s(z^,z^)  - gCz^.z^)]  = 


(z  - z ) g(z  ,z  ) 
1 2 ^ 1’  2 


(A.  7) 
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* 


one  finds 


X 


.5=o)exp  iq  (z 
1 ^ z 1 


(A.  8) 


Ignoring  the  terra  ^ d( | Vq( z)  | ,0)dz  in  the  exponential  which  can  only 
make  the  integral  sraaller,  changing  to  variables  u = z^-  z^  and 
i Zg),  approximating  their  limits  by  -»  < u < » and  0 < T < R 

and  setting  g{z^,  z^)  « g(z,z}  then  yields 


<I>  - <i>^ 


Rp(O.O) 

3L^ 


2 

R<y  > 
3L 


(A. 9) 


The  correction  to  Markov  approximation  to  <(^  (2)  ^il)>  can  be 
analyzed  in  the  same  way.  It  is  fractionally  small  as  long  as  < y^>  R/L 
is  small. 


f 


I 
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APPENDIX  B CORRECTIONS  TO  GAUSSIAN  STATISTICS 


/ 

To  begin  with  something  simple,  consider  (1  It  is  given  by 
the  path  integral  in  Eq.  (3.1)  with  t = t and  the  end  point  conditions 

-»  -♦  4 ^ 

r (0)  = r (R)  = 0.  Changing  variables  to  u , u , v and  v defined  by 
i i 1.  2 X « 


rj^(z)  = u^(z)  + iv^(z)  + Ju^Cz)  + 

T^(z)  = u^(z)  + iv^Cz)  - “ iv^(z) 

r (z)  = u (z)  - iv  (z)  + iu  (z)  - Jv  (z) 

3 X 2 2 X 

r^(z)  = u^(z)  - iv^(z)  - * iv^(z)  (B.l) 

the  integral  over  u (z)  can  be  done  and  it  produces  a 6 -functional 

which  forces  u (z)  to  vanish  everywhere.  With  the  end  point  conditions 
-»  -♦  -* 

u (0)  = u (R)  = 0,  the  only  solution  is  u (z)  = 0.  The  quadruple  path 
2 2 2 

integral  then  reduces  to  the  double  path  integral  over  paths  v and  v 

X 2 


K 

^ (paths)exp[ik,yv'(z)  • v'(z)dz  - 

1+k  *■  0 ^ ^ 


(B.2) 


where  M was  defined  in  Eq,  (3.2)  and  for  u^  = 0 it  is  explicitly 

Ri 

M = y l2d(lv^(z)l  ,o)  + 2d(lv2(z)l,0  )-  d(lv^(z)  + V2(z)l,o) 
0 

- d|lv^(z)  -V2(z)l,o)  dz 


(B.3) 


The  two  regions  (a)  and  (b)  discussed  in  Section  III  are  iVj^(z)|  < L/4 
“>  !-♦  1 

with  V arbitrary  and  jv  (z)j  < L/^  with  v arbitrary. 

2 2 X 

It  was  pointed  out  in  the  text  that  M is  of  order  unity  or  smaller 
throughout  regions  (a)  and  (b).  Actually,  there  is  a further  region 
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(having  some  overlap  with  (a)  and  (b))  where  M can  be  small.  It  is  (c) 
Iv  (z)|  < L//9  , |v  (z)|  < L//5  and  owes  its  existence  to  the  fact  that 

X dt 

■ ->  1 I 1 ^ 

when  both  |v^(z)|  and  jv^Cz)!  are  small,  M is  quartic  in  the  v's.  In 
- -M 

all  other  regions  of  path  space,  e is  exponentially  small. 


Our  first  task  is  to  dispose  of  the  extra  region  (c)  by  showing 
that  for  small  o the  volume  of  path  space  occupied  by  this  region  is 
exponentially  small  compared  to  the  volumes  occupied  by  regions  (a) 
and  (b).  An  estimate  of  the  volume  of  path  space  occupied  by  region  (a) 
is 


— f(^ {'pa.ths')  exp  [ik^yv'(z)  • v'(z)dz (z)) 

kt  L 0 ^ 2 L R 0 1 


dz 


= 1 (B.4) 


where  the  integral  is  done  by  integrating  over  v^Cz)  which  produces  a 
6-functional  that  forces  v^(z)  to  vanish.  An  estimate  of  the  volume 


occupied  by  region  (c)  is 


1 /*  2 [ A, 

(paths)  exp|ik_/v^ 


0 

(6/a) 


2 (3  2/3 

sinh  ■*/—  + sin  *)  — 

Jfa  la 


(z)  • v^(z)dz 

2KL 


-W- 

24  Ifa 
— e 
O' 


^V^(z)  + V^Cz)) 


dz 


(B.5) 


where  the  value  of  the  path  integral  is  taken  from  Ref.  (6).  For  small 

a,  the  volume  occupied  by  region  (c)  is  therefore  exponentially  small 

compared  to  the  volume  occupied  by  regions  (a)  and  (b).  This  result, 

which  may  surprise  some  readers,  deserves  an  explanation.  In  region  (a) 

where  jv  1 is  always  less  than  L/$ , the  factor  exp  mt  fv'(z)  • v*(z)  dzl 
1 1 2 J 


= exp 


(z)  • v^^(z)dz 


will  restrict  I vj 


\ to  values  less  than 
# /(kLR) . Typical  values  of  jv  | will  then  be  R^|v^^l/6  L/of.  At  a 


Ik 


given  range  point  z , the  variables  v (z  ) and  v (z  ) span  a four- 
o 1 o 2 o 

dimensional  space.  In  this  space  the  volume  occupied  by  paths  in 

4 2 

region  (a)  is  roughly  L /(ia)  , At  the  same  point  the  volume  occupied 
.V  4,-2 

by  paths  in  region  (c;  is  roughly  L /*  . Thus  at  each  point  z , the 

2 ° 

volume  associated  with  region  (c)  is  a factor  of  a smaller  than  that 

occupied  by  region  (a).  To  compute  the  total  volume  in  path  space, 

one  has  to  multiply  together  the  volumes  at  each  range  point  z , taking 

o 

into  account  the  fact  that  the  paths  cannot  bend  too  rapidly.  The  path 
integrals  in  Equations  (B.4)  and  (B.5)  do  just  this.  The  resulting 

exponential  ratio  of  volumes  should  no  longer  be  a surprise  since  at 

^ 2 
each  range  point  the  ratio  is  down  by  a . 

It  is  therefore  sufficient  to  consider  only  paths  lying  in  regions 

(a)  and  (b) . The  fact  that  integrating  separately  over  regions  (a)  and 

(b)  leads  to  a slight  over-counting  can  also  be  Ignored.  The  volume  in 

path  space  where  regions  (a)  and  (b)  overlap  is  even  smaller  than  the 

volume  occupied  by  region  (c) . Now  as  was  pointed  out  in  the  Section  III, 

(a) 

for  most  paths  in  region  (a)  M M where 

o 


= 2^(lv^(z)l,o)( 


(B.6) 


and  for  most  paths  in  region  (b)  M « M with 

o 


K 

”o*^^  = 2 yd[lv2(z)l,o)ds 


(B.7) 


The  path  integral  in  Eq.  (B.2)  can  then  be  expanded  according  to 
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R 


(paths )exp 


(z)  • \*^{z)dz  - M 


, V (a) 

(z)dz  - M 
2 o 


/ (a)  im 
(“o  -") 


CO 

^ ^ y*d^(paths)expfily  v'(z)  • 
m = 0 1-0 

oo  p / (b)  \in 

(paths  )exp  ik|v^(z)  • v'(z)dz  - M — - — ; (B 

_ I / 1 2 O D1  • 

m = 0 I •& 


8) 


which  is  an  asymtotic  series  in  or.  It  is  not  a convergent  series  because 
(exponentially  small)  contributions  from  region  (c)  and  the  overlap  of 
regions  (a)  and  (b)  are  not  being  treated  correctly.  The  m = 0 terms 
correspond  to  Rayleigh  statistics  and  the  m = 1 terms  are  the  first 
correction.  They  will  be  computed  explicitly  below. 


Eq . (B.8)  generalizes  to  an  arbitrary  correlation  in  the 

obvious  way.  For  a 2nth  order  moment  there  are  n ',  important  regions 

of  path  space.  In  each  such  region  there  is  an  M given  by  the  analog 

o 

of  Eq  . (3.^)  or  (3* 5)-  The  generalization  of  Eq.  (B.8)  is  then  a 

sum  of  n'.  terms,  each  of  which  is  a series  of  powers  of  the  appropriate 


M - M. 
o 


The  path  integrals  for  the  m = 1 terms  in  Eq.  (B.8) 

ated  by  inserting  a spectral  representation  for  M - M, 

o 

three-dimensional  Fourier  transfer  of  p (See  Eq.  (5.3)), 


can  be  evalu- 
If  p is  the 
then  in  region 


(a) 


cos 


(B.9) 
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and 


fl 

^^^d^Cpaths)  exp|ik  J‘\ '^{z)  • {z)  Az  - ^^^^j  - 

n O O 


M = 


R r ® 

^ Jdz'J'd  q p (|  q|  ) J' d (paths)  exp  j” ’ ( ''i ' ~ k ^ 


dz 


X 1 - cos  ^ q • (z  ^)) 


(B.IO) 


Since  M depends  only  on  v , the  integral  over  v can  be  done  and  it 
o 1 ^ 

produces  a 6-functional  which  forces  v' ' (z)  - ^ 6(z  - z")  to  vanish. 

1 ^ 

In  terms  of  the  Greens  function  g(z,z')  defined  by 


— - g(z,z')  = 6 (z  - z') 
3z 


(B.ll) 


with  boundary  conditions  g(0,z')  = g(R,z')  = 0,  v^Cz)  is  constrained 
to  be  ^ g(z,z').  In  ^ and  cos  ( q * v (z^)|  one  can  then  set  v (z) 

q , \ 1 / 1 

equal  to  - g(z,z  ) and  the  remaining  path  just  gives  (l) . The  calcula- 
tion of  the  correction  in  region  (b)  is  identical  and  to  leading  order 
in  a 


<I  > - 2<I)"  ^ 
<.)^  ■ 


4nk 


q p(| ql)  Q(z,  Iqj) 


(B. 12) 


with 


Q(2,1^  ) = 2 


1 - cos 


R 

( 1“  g(z.  z))  exp|-2  g(z,z')  ,oj  dz' 


(B.13) 
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An  examination  of  the  integral  jon  the  right-hand  side  of  Eq.  (B. 12) 
shows  that  for  small  or,  (i)  d g(z,  z ^) , o)  can  be  approximated  by 

, (ii)  2 1 - cos(^g(z,z))  can  be  set  equal  to 

y— g(z,z)y  and  (iii)  the  dominant  contribution  comes  from  the  regions 

Z 0 and  z « R where  g(z,z^)  » z(R  - z^)/R  and  z*(R  - z)/R.  The  con- 
tributions from  the  regions  z » 0 and  z R are  the  same  and 


= 8 Tik^y dz  J'd^q  P ( 1 q 1 )( 


exp  - 


2-*2 

Rp(0,0)z  q 


- 1‘ 


L Jq  P(q,0)dq 



00 

yqp(q,0)dq 

0 


(B.14) 


The  correction  to  (!*')  involves  nl  regions  of  path  space  and  in 

each  of  these  regions  there  are  n(n  - l)/2  terms  in  M - M which  differ 

o 

only  by  permutations  of  the  paths.  The  result  is  that  the  correction  to 

IV  2 

(l  > is  n;(n)(n  - l)/4  times  the  correction  to  (l  ). 

Moving  on  to  a more  complicated  object,  consider  (l(l)T(2)).  It  is 

given  by  the  path  integral  in  Eq.  (B.2)  but  the  end  point  conditions  on 

the  paths  are  now  v, (0)  = v (R)  = 0 and  v^(0)  = r - r , v (R)  = r - r 
1 1 2 olo2  2 12 

and  now  M is 
R 

M =^|2d  (1  v^(z)  1 ,o)  + 2d^lv2(z)|  ,t^  - t^  - d (1  v^(z)  + v^Cz)!  ,t^  - t^) 
0 ^ 

- d (l v^(z)  - v^Cz) I , - t^)  dz  (B. 15) 


The  integration  over  region  (a)  gives  (l)  (1  + corrections)  while  the 

2 -D(l  2) 

Integration  over  region  (b)  gives  (l)  e ’ (1  + corrections).  As 
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Indicated  the  corrections  in  region  (b)  are  proportional  to  e and 

are  a small  effect  of  no  particular  consequence.  The  corrections  in 

-D 

region  (a),  on  the  other  hand,  are  small  but  do  not  contain  e and 

hence  fall  much  less  rapidly.  This  leads  to  a coherence  tail  in 

(l(l)I(2))  which  is  not  present  in  , The  interesting 

corrections  in  region  (a)  are  computed  by  changing  variables  from  v (z) 

and  v„(z)  to  v (z)  and  w(z)  = v (z)  - - (r  - r ) — (r  , " r „) 

2 1 2R12RoXo2 

and  then  proceeding  in  exactly  the  same  way  as  before.  The  result  is 

<I(1)I(2)>  = <I>^[1  + e + Y(r  - r ,r  “ r t,  " t,)] 

1 2 ol  o2  1 2 

(B.16) 


where 


Y(r.r  , t)  = 2Trk 
o 


2 

d q exp 


zq  • r (R  - z)q  • Tc 

i + i ^ 


i:(lql  , t)Q(z,  jql) 
(B.17) 


and  Eq.  (B.12)  is  not  reproduced  for  1(1)  = 1(2)  because  a small  term 
-D 

of  order  ye  from  region  (b)  has  been  dropped.  For  small  a this  integral 
can  be  simplified  in  the  same  way  that  Eq.  (B.14)  was  obtained  from 
Eq.  (B.12).  It  becomes 


Y(r,r  ,t) 
o 


00 

L /q^P(q,t)[j  (qlrj)  + J (qjr  j ) ]dq 
./  o o o 

_Q 


( 

/ 


qp(q,0)dq 


(B.18) 


0 


Corrections  to  more  complicated  correlations  and  terms  of  order 

2 

or  or  higher  can  also  bo  computed  — the  only  obstacle  being  the  labor 
involved.  The  calculation  of  the  general  coherence  tail  involves  only 


some  combinatorics.  It  is 


n 

- I 

K,J=1 


{B.19) 


where  Y(k  - j)  = y(?j^  - ) , the  "Gaussian  terms"  are 

what  one  would  compute  from  the  Gaussian  distribution  and  all  terms  of  order 
e ^ have  been  dropped. 
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APPENDIX  C;  CORRECTIONS  TO  GAUSSIAN  STATISTICS  FOR  p < 2 


Eqs.  (B.12),  (B.13),  (B.16)  and  (B.17)  of  App.  (B)  do  not 
assume  a single  scale  media  and  will  be  the  starting  point.  For 
p < 2,  in  either  the  fully  or  the  partially  saturated  regime  Q 
can  be  approximated  by 


Q(z, |q| ) 


g(z,2) 


2 


exp 


ip 

IqI  g(z,z) 

P 

P+1 

kL 

(C.l) 


where  the  cosine  has  been  expanded,  the  short  distance  expansion 
for  0 Eq.(7.2)  has  been  used  and  the  identity 


R 

I |g(z,z' ) l^dz’  = ^^lg(z,z)lP 

0 


(C.2) 


has  been  employed. 

In  the  fully  saturated  regime  the  main  contribution  to 
Eqs.  (B.12)  and  (B.17)  again  comes  from  z ' 0 and  z ““  R.  Using 
Eq.  (C.2)  for  Q then  yields  Eqs.  (B.14)  and  (B.18)  with  a re- 
placed by  a'  where  a'  is  defined  in  Eq . (7.3). 

In  the  partially  saturated  regime  all  values  of  z 
contribute  to  the  integral  but  the  dominant  contribution  comes 
from  large  |q|  where 

g(q,t)  = g(0.t)2P  (r(l+^p))^  _sin(Trp/2) 

4it^  iLlPlqlP"^^ 


For  n = 2,  Eq.(8.3)  with 

C(p)  = 2P  (r(UAp)r(p -1))^  sin(TTp/2)r  ((4  - p)/p) 

PTT  6^^~P^  r(2p-2)  (C.4) 
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is  obtained  by  inserting  Eqs.(C.3)  (with  t=0)  and  (C.l)  in 
Eq.(B.12).  The  extension  to  general  n works  in  the  same  way 
as  before.  The  coherence  tails  in  the  partially  saturated 
regime  are  obtained  by  inserting  Eqs.(C.3)  and  (C.l)  into 
Eq,(B.17).  For  r,  ^ 0 this  leads  to  integrals  which 
cannot  be  done  analytically. 
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APPENDIX  D:  THE  PARTIALLY  SATURATED  REGIME  FOR  p > 2 


It  is  difficult  to  make  quantitative  statements  about 
the  partially  saturated  regime  when  4 > p > 2.  There  is  how- 
ever some  qualitative  information. 

Eq.(1.14)  holds  and  <<f  > is  equal  to  ^exp  in 

all  regimes  as  long  as  the  Markov  approximation  is  valid. 

Furthermore  (w)^  continues  to  be  given  by  Eqs.(1.19) 

- (1.21).  The  argument  that  any  correlation  involving  an  un- 
equal  number  of  ^'s  and  S 's  vanishes  also  goes  through  as 
before.  Thus  ^ is  uniformly  distributed  in  phase.  The 
difficulty  arises  when  one  attempts  to  compute  the  non-vanish- 
ing higher  moments.  ' 

The  statistics  are  not  Gaussian.  This  can  be  verified 
by  assuming  that  they  are  and  then  computing  the  corrections. 

They  are  not  small.  Some  information  can  be  obtained  however  j 

I 

by  comparing  the  path  integral  for  ^ (2))^  (<S^(1))^^  with  ^ 

that  for  ^ I ^(2)1^  |(^(1)|^)  • Upon  doing  this  one  finds  that  j 

always  phase  wraps  as  shown  in  Fig. (5b).  It  turns  out  that  i 

the  typical  space  time  scales  over  which  the  phase  and  intensity 
change  are  those  listed  in  Table  D-1  (the  parameter  6 was  de- 
fined in  Eq . ( 8 . 4 ) ) . 


* 
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Phase 

L/'t 

T/4> 

£n  I 

Ez2 

2 

(T/$)(n/$)^‘P 

(T/<1.)(R/«>)'*“P 

p - 6 > 2 

2 > p - 6 > 0 

Scale  Length 

Scale  Time 

Table  D-1.  Time  and  Space  Scales  Associated  with 

Phase  and  £n  I in  the  Partially  Saturated 
Regime  with  4 > p >2. 


4 / p 

Note  that  for  partial  saturation  where  •t  > 1 but 

<I>/n  < 1 the  rate  at  which  the  intensity  changes  is  always 
small  compared  to  the  rate  at  which  the  phase  changes.  Ex- 
amining more  complicated  correlations  leads  to  the  conclusion 
that  at  a fixed  frequency  can  be  represented  as 

(^(j)  = SqU)  exp[i  4>(j)]x(j) 

where  (})(j)  is  a real  Gaussian  random  variable  with  ^ 4>(3))  ^ 
and 

((^(i)  - 4i(j))^)  = D(i,j) 

The  other  factor  x is  an  independent  (of  (}>  ) complex  random 
variable  about  which  only  three  things  are  known: 

(1)  any  correlation  involving  an  unequal  number  of  x's  and 

♦ 

X 's  vanishes 

(2)  <lx|^>  = 1 

and 


(D.l) 

= 0 

(D.2) 
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(3)  the  decorrelation  lengths  and  times  for  x those  listed 

under  intensity  in  Table  D-1. 

To  see  what  the  representation  In  Eq.  (D.l)  means  consider 

( x*(2)  x(i)) 


= exp  [-iD(l,2)]  < |x|^ 


(D.3) 


where  to  get  the  second  line  one  notes  that  ^X*(2)  X(l)) 
will  be  approximately^  1x1^^  for  all  space  or  time  separa- 
tions such  that  exp  [ “§Dj  is  not  vanishingly  small.  Thus, 
Eq.(1.19)  is  reproduced,  as  it  should  be.  Similarly, 


(1))^) 


^ [-20(1,2)]  <1x1^  ) 

(f^(2)r  (6  o(i)r  j \ / 


(D.4) 


and  this  correlation  is  known  up  to  a constant.  However, 
all  that  is  known  about  the  intensity  correlation 


<I(2)  Kd)  = ( |x(2)|2  Ix(l)l^) 


(D.5) 


are  its  space  and  time  scales, 


The  extension  to  unequal  frequencies  is  straightforward. 

(I  \ 2 

----  - 1 

u)  j and 

* g ' 

(<^' ) - m').  The  higher  order  moments  of  x 

are  again  non-Gaussian  and  unknown.  However,  their  width  in  m 

is  large  compared  to  gj  . As  in  Sec.(8''  this  means  that  pulse 

O 

propagation  is  dominated  by  wander  rather  than  spreading. 
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As  was  mentioned  in  Sec. (8),  there  is  a case  where  the 
non-Gaussian  statistics  of  x can  be  studied  in  detail.  It  is 
for  correlations  in  time  when  p = 2 and  6=0  and  is  explained 
in  Ref . (7). 

The  above  results  are  most  easily  derived  using  the 
Fermat  path  formalism  of  Sec. 6.  One  can  work  out  the  joint 
probability  that  two  paths  will  satisfy  the  perturbed  ray  equa- 
tion. In  the  partially  saturated  regime  with  p > 2 it  turns  out 
that  the  Fgrmat  paths  are  highly  correlated  and  tend  to  lie  with- 
in of  each  other.  Studying  averages  of  y and  y'  along 

correlated  Fermat  paths  then  leads  to  the  above  conclusions.  The 
detailed  calculations  are  relatively  straightforward  but  tedious 
and  will  not  be  given  here. 
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APPENDIX  E;  CORRECTIONS  TO  THE  MARKOV  APPROXIMATION 
FOR  INHOMOGENEOUS  ANISOTROPIC  MEDIA 

If  the  X dependence  of  p is  evaluated  along  the  unper- 
turbed ray  then  the  first  correction  to  the  Markov  approximation 
can  be  evaluated  for  a general  homogeneous  anisotropic  medium. 
Let 


p(x,  t;  s(z)  + e z)  = I d^Jle^^  ^p(l,t;z) 
z 


then  the  generalization  of  Eq.  (A. 6)  is 


(E.l) 


R R 

= 2k^  I dz^ 

0 0 


dZ2  I d^q  dq^p[q  +l^q^,0;  i 


exp 


X sin 


i(q2+  1 ’S'  (i(Zj^+  z^))  (z^-  z^) 


R 


sin 


iq-VQ(Zj^) 


iq-Vo(z2)jexp 


jexp  - I d(vQ(z) ,0;z)dzj  (E.2) 


where  q = (qj^.qy)  is  a two  dimensional  vector,  d is  defined 
in  Eq . (9 . 11 ) and 


'o^^>i  = 1?  - g.jCz.Zg)! 


(E.3) 


with  g^j  defined  in  Eq.(9.22). 

For  an  isotropic  medium  Eq.(E.2)  can  be  analyzed  in  the 

same  way  as  Eq.(A.6)  and  one  finds  that  (<I>  - <I>)/  <I> 

0 0 

is  of  order  of  the  r.m.s.  multiple  scattering  angle. 

It  is  also  straightforward  to  analyze  Eq.(E.l)  for  a homo- 
geneous but  anisotropic  medium.  The  result  of  doing  this  was 
stated  in  Sec.(9A). 
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APPENDIX  F;  CORRECTIONS  TO  GAUSSIAN  STATISTICS  FOR 
INHOMOGENEOUS  ANISOTROPIC  MEDIA 


When  the  approximation  of  Eq.(9.7)  for  the  correlation 
between  two  paths  is  made,  it  is  possible  to  compute  the  cor- 
rections to  Gaussian  statistics  in  the  saturated  regimes.  The 
calculation  is  a fairly  straightforward  generalization  of  that 
done  in  App.B  and  only  the  final  result  will  be  given. 

Define  a function  p,(q;z)  where  q = (q  ,q  ) by  P,(q,t:z)  = 

X X y j. 

p(q  - e (s ' (z)  "q ) , t ; z)  where  p is  defined  in  Eq.(E.l).  Then  the 
z 

analog  of  Eq.(B.12)  is 


d q p^(q,0;z)  Q^(z,q) 


(F.l) 


and  d and  g. . are  defined  in  Eqs.(9.11)  and  (9.22). 

Using  D and  <J>  as  defined  by  Eqs.(9.21)  and  (9.6)  it  is 
possible  to  show  that  the  right-hand  side  of  Eq.(F.l)  is  small 
in  the  fully  saturated  regime  and  in  the  partially  saturated 
regime  for  p < 2.  As  before,  in  the  fully  saturated  regime  the 
dominant  contribution  comes  from  the  regions  z ?:  0 and  z * R 
and  in  the  partially  saturated  regime  p^  can  be  approximated 
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by  its  asymptotic  form  for  large  |q|.  Also,  can  be  simplified 
by  expanding  the  cosine  and  replacing  d by  its  expansion  for 
small  The  detailed  calculation  which  is  then  fairly  straight- 

forward will  be  left  to  the  reader. 

The  generalization  to  <l”>  works  in  the  same  way  as  in 
App. (B) . 

The  coherence  tail  is  given  by  Eq.(B.16)  with 

^01-^02’ ^ ^ 

R 

= 2Tik^  dz  d^q  exp  iq-t!kz)  p^^Cq , t ; z)Qj^(q , z)  (F.4) 

where  ^ is  defined  in  Eg. (9. 14).  The  aporoximat ions  mentioned 
above  can  also  be  made  in  the  integral  for  y • Finallv.  Eq.(B.19) 
holds  with  Y given  by  Eq.(F,4). 
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FIGURE  CAPTIONS 


I 


FIG.l  The  difference  between  multiple  large  angle  scattering 

(a)  and  multiple  small  angle  scattering  (b). 

FIG.  2 Parameter  regimes  in  <I>  - space. 

FIG. 3 Illustrating  Eq.(1.18).  The  signal  n ' will  lie,  with 

90%  probability,  within  the  circles;  (1)  for  $t/T 
small,  (2)  for  iJ>t/T  ~ 1 and  (3)  for  $t/T  large.  The 
location  of  the  signal  n at  t=0  was  an  unlikely  one 
lying  outside  the  90%  probability  circle  for  a 
Rayleigh  distribution. 

FIG. 4 A path  in  the  path  integral  for  n=6. 

FIG. 5 (a)  The  schematic  track  of  a signal  satisfying  Gaussian 

statistics  in  time. 

(b)  The  track  of  a signal  which  moves  faster  in  phase 
than  amplitude  (phase  wrapping). 

FIG. 6 Propagation  through  an  anisotropic  medium.  The  blobs 

are  schematic  inhomogeneities  and  the  heavy  directed 
lines  are  the  unperturbed  propagation  path  at  various 
angles  with  respect  to  the  long  axis  of  the  blobs. 

FIG. 7*  Propagation  in  a channel.  The  channel  axis  (z-axis) 
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